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ABS'ERACT 


Dyadic correlation properties of a class of second degree 
binary Boolean functions, called NON-REPETITIVE QUADRATIC FORMS 
(NRQFs) / are analysed using Walsh-Hadamard transform techniques. 
It is shown that these 2^-point functions of n variables possess 
ideal dyadic auto- correlation and asymptotically ideal cross- 
correlation properties. Walsh-Hadamard spectral coefficients 
of NRQFs are found to be of uniform magnitude. It is proved 
that all 2^-point functions with uniform magnitude spectral 
property possess ideal, dyadic auto- correlation and vice-versa. 
Such functions are shown to be optimal as perturbation signals 
for identification of linear dyadic-shift invariant systems. 

The class of these functions, known as bent functions, are 
presented briefly with their properties. The NRQFs constitute 
a sub-class of binary bent functions. Some interesting results 
on dyadic correlation properties of two-dimensional arrays, 
synthesised using NRQFs and Walsh functions, are derived. For 
comparison with linear- feedback shift-register (DFSR) sequences 
which are known to possess good periodic and aperiodic correla- 
tion properties, periodic and periodic correlation parameters 
of some NRQFs are computed. Minimal LFSR circuit realiziatLons 
of NRQFs are studied and their complexity ccanpared- with maximal- 
length DFSR sequence generators^ Finally# a method is suggested 
for design of binary segumces wu-tti desired dyadic correlation 



CHAPTER I 


INTRODUCTION 

Considerable research has gone into the analysis and 
synthesis of sequences with good correlation properties 
[soehmer, 1967; MacWilliams and Sloane, 1976; Shedd and 
Sarwate^ 1979; Sarwate and Parsley/ 1980 ]. The resulting 
achievements have been both a mathematician's delight and an 
engineer's boon. Binary and poly-phase sequences/ both have 
attracted researchers. However# the focus of attention have 
been binary sequences as they are easily generated using the 
available hardware. Two level sequences with good cyclic ^d 
linear correlation properties find many applications in their 
respective domains. From the more abstract systems' viewpoint# 
correlation# defined in a given way# is closely related to 
systems characterised by input-output convolution relations 
defined in a similar way. Further# correlation analysis can 
be performed in the transform domain using the system- specific 
transform. Consequently, both the sarrple domain and the trans- 
form domain properties of sequences are of use in their corre- 
lation studies. 

Binary sequences synthesised using linear- feedback shift- 
register (LFSR) configxirations are mathematically elegant# 
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easily generated and rich in their potential for use in various 
applications [ Golomb/ 1964; Golomb, 1967] . A sub-class of 
these, called maximal-length sequences (m-sequences) or pseudo- 
random binary sequences (PRBSs) , possess asymptotically ideal 
linear (aperiodic) and cyclic (periodic) correlation properties. 
They have found extensive use, traditionally, in Radar Ranging 
[Golomb, 1964j , Cryptography. [Geffe, 1973; Feistal, et al,,l975j 
and Scrambling [Henriksson, 1972], and, more recently, as 
wide-band signature sequences in Spread-Spectriam Code- Division 
Multiple-Access (SS-CEMA) communication systems [ Pursley and 
Sarwate, 1977; Pursley and Roefs, 1979] . Actually, large sets 
of sequences, known as Gold [Gold, 1968 ] and Kasami [ Sarwate 
and Pursley, 198o] sequences derived from a combination of 
m~sequences are used for SS-CIMA. systems, nv-sequences, owing 
to their good auto- correlation property, constitute; a near- 
ideal set of signals for identification of linear cyclic-shift 
invariant (LCSIV) systems. Cyclic convolution and discrete 
Fourier transform (DFT) relate the input and output of these 
systems in the sample and transform domains, respectively. 

DFT row vectors are their ei gensi goals. Hence, cyclic correla- 
tion analysis of LFSR and other sequences can be performed 
using DFT techniques. Non-binary poly-phase sequences possess- 
ing ideal cyclic correlation are also known. These are the 
Frank-Zadof f-Chu, (FJSC) sequ^ces [ Frank and Zadoff, 1962; Chu, 
197 ^, 
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Some of the well known sequences having near-optimal 
aperiodic correlation properties are Barker (binary) and 
Huffman (non-binary) sequences ^uffman^ 196^. Unfortunately/ 
binary Barker sequences are known to exist only up to length 13. 


Generalized Barker sequences |j3olomb and Scholtz, 

Huffman sequences are of not much use practically. Hence, 
other sub-optimal binary sequences - viz. pulse compression 


r 


codes ,|Boehmer, ^ even- shift orthogonal sequences j^aki 

and Miyakawa, i960 etc. -have also been studied. 


1. 1 SCOPE OF THE PRESENT WORK 

The present work is devoted to correlation analysis of 
a class of binary sequences with good dyadic correlation 
properties. Dyadic correlation analysis uses Walsh-Hadamard 
transform (WHT) techniques because Walsh functions are eigen- 
signals for linear dyadic- shift invariant (LDSIV) systems 
whose output is dyadic convolution/correlation between the 
input and the system sanple vector. Accordingly, WHT techniques 
are applied to a class of n variable 2 '^-point binary Boolean 
functions called NON-REPETITIVE QUADRATIC FORMS (NRQFs) 
^arpovsky, 1970 / for study of their dyadic correlation 
properties. WHT coefficients of NRQFs are computed by 
correlating them (NRQFs) with linear Boolean functions (DBFs) , 
as each row of an N xN (N = 2^) Walsh-Hadamard matrix is one 
of 2^ LBFs of n variables. The WT coefficients turn out to 
be of uniform magnitude, 2^^^. Uniform magnitude sequency- 
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spectriam results In the 2^-point |^+1» -valued NRQFs having 
ideal dyadic auto- correlation. Thus/ they constitute an ideal 
set of perturbation signals for identification of LDSIV systems. 
NRQFs constitute a sub-class of another known set of functions, 
called BENT FUNCTIONS/ which have the uniform spectral property. 
These have also been presented briefly. 

Further, cross-correlation properties of NRQFs are 
derived. They are shown to be asymptotically ideal. This 
has been achieved using the permutation cycle structure 
relationships between NRQFs. An important result relating 
cross-correlation function of two NRQFs to the cross-correlation 
functions of each of these with a third NRQF, is derived. 


In order to compare the correlation performance of NRQFs 
vis-a-vis m- sequences, periodic and periodic correlation 
parameters of some typical NRQFs are computed. However, they 
turn out to be considerably worse. Another parameter for 
comparison is the complexity of minimal LFSR circuits as NRQF 
generators vis-a-vis m-sequence LFSR synthesizers. Number of 
stages of minimal LFSR circuits needed to generate NRQFs of 


2/ 4, 6 and 8 variables are conputed using Massey's shift- 
register synthesis algorithm |Massey, I960. 


Taking a cue from the one- to-one relationship between 
existence of certain integer difference sets and sequences 
with good cyclic auto— correlations (nv-sequences etc.) 

1969 1, a method is suggested for design of sequences with 




5 


desired dyadic auto- cor relation properties using what have 
been defined as dyadic difference sets. 

Finally# dyadic correlation properties of NRQFs and Walsh 
functions are utilized for synthesis of two-dimensional 
hi, - 1 } -valued arrays with good correlation properties. Two 
methods of synthesis are presented and the resulting correla- 
tions derived. A judicious choice of NRQFs and Walsh functions 
as the component arrays leads to som.e interesting results. 


1. 2 ORGANISATION OF THE THESIS 

In Chapter 2 , we present an integrated brief review on 
sequences with good periodic and aperiodic correlation proper- 
ties. Beginning with the definition of periodic correlation# 
we relate it to the DFTs of the component arrays and derive 
certain conditions to be satisfied by their DFTs for the arrays 
to possess ideal auto- and cross- correlations. FZC and Nth 
roots of unity - sequences are presented next as specific 
instances of sequences with good periodic correlation properties- 
Then, we go on to a brief exposition of m- sequences - their 
generation and some iirportant properties. Their as ynp to tic ally 
ideal auto- correlation is related# next, to their DFTs. An 
overview of irrportant properties of LGSIV systems and a scheme 
for their identification are given. The chapter ends with a 
brief review of sequences with good aperiodic correlation 
properties. 



In Chapter 3, we present an overview of dyadic correla- 
tion and LDSIV systems. Dyadic correlation/convolution and 
their relation to LDSIV systems is given which is followed by 
some of their elementary properties in terms of the Walsh- 
Hadamard transform. A scheme for the identification of LDSIV 
systems is presented next. A typical dyadic correlator 
configuration used in the scheme is also presented. Finally, 
a note on WHT-based classification of Boolean functions is 
added as an aid to dyadic correlation analysis in Chapter 4. 

Chapter 4 is a detailed presentation of non-repetitive 
quadratic forms (NRQFs) . Beginning with their definitj.on, 
elementary properties and method of generation using counters, 
minimal LFSR circuit realisations of some typical NRQFs are 
tabulated next. Then, WHT coefficients of NRQFs are derived. 
These are used to analyse their auto- and cross-correlation 
properties in detail. Finally, aperiodic and periodic correla- 
tion parameters of a few NRQFs of 4, 6 and 8 variables are 
presented. 

Bent functions, a generalisation of NRQFs, are reviewed 
in Chapter 5, Inportant properties of these functions are 
presented. Scope for enumeration of bent functions of a given 
number of variables is also analysed. Quadratic bent functions 
are given as a special class of bent functions which can be 
cDirpletely and easily enumerated. The chapter ends with 
proving the ideal auto- correlation property of bent functions. 



7 


Chapter 6, the concluding chapter, summarises the 
achievements of this work. Certain related directions for 
further endeavours are suggested. Further groundwork is 
needed for these to lead to certain conclusive results. 



CHAPTER 2 


SEQUENCES WITH GOOD PERIODIC AND 
APERIODIC CORRELATION PROPERTIES 


This chapter is a review of linear cyclic- shift invariant 
systems/ DFT and sequences with good periodic and aperiodic 
correlation properties. Beginning with a definition of cyclic 
correlation in Section 2.1 and introduction of the DFT approach 
to its analysis/ conditions on DFTs of sequences with ideal 
correlations are derived in Section 2.1.1. Section 2*1.2 gives 
FZC and Nth roots of unity - sequences as instances of sequences 
possessing the required DFTs for ideal auto-correlation and 
ideal cross-correlation, respectively. Section 2.2 is a brief 
exposition of m-sequences - their generation and a few important 
properties. Their asynptotically ideal auto-correlation and 
the required DFTs are also presented. This DFT approach to 
m-sequences is rare in literature, LCSIV systems, their 
elementary properties and a scheme for tbeir identification 
using m-sequences constitute Section 2.3. Finally, Section 2.4 
with a brief overview of sequences with good aperiodic correla- 
tion properties marks the end of this chapter. 
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2.1 CYLIC/PERIODIC CORRELATION 

We begin with the preliminaries of defining cyclic or 
periodic correlations. In general, given two N-length complex- 
valned sequences: 


... 


“n-i^ 


2 = 




the cyclic shift by s of u is 


T®u = (u„ u 


N-s' N-s+1 


"^O' .^1' * * •'^N-l-s^ 


i.e. (T^u) . = u. , (i-s) mod N. 

***** i JL S 

The periodic auto- correlation for any cyclic shift is defined as; 

N-1 

Pn =^u, T®u E u. , (i-s) mod N, 

u ^ i=0 ^ 


0, 1, .. . (N-l) , 


( 2 . 1 ) 


where ^ denotes the inner product of the conponent sequences 
and is the comp lex- conjugation. 

Similarly, cyclic cross-correlation between u and y is: 

N-l 

^uv^®^ ® ^ ^i'^i-s ' 

i^O 


= 0, 1,... (N-l) 


( 2 . 2 ) 


As is evident from the definitionj 


P^(s +N) = P^Cs) and p^^(s +N) = P^^(s) 


(2.3) 
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Of prime interest here are ideal and asymptotically ideal 
(near-ideal) correlations and the sequences possessing them. 

A sequence u is said to possess ideal periodic auto- correlation 
if; 


P^(s) 



(2.4) 


i.e. if the sidelobes are identically equal to zero. 


Ideal cross- correlation between u and v inplies that: 


p^^(s) = 0 , for all s = 0, 1, (N-l). (2.5) 


The requirements for these ideal correlations can also 
be translated into certain conditions to be satisfied by the 
DFTs of the coinponent sequences. The DFT of an N-length array 


X = (x^x^ 


^-1 


) 


T 


is defined as: 


X(s) 



N-l 

E 

i=0 


% 

Xj^ (- j 


2 ^ si . 

N I' 


s =0, 1,...,(N-1) 


In the matrix form# this is: 



1 

f 


F 


( 2 . 6 ) 


(2.7) 



1 X 


where 


L 





a 

2 

a 

N- 

2 

a 

4 

a 

N- 

N-1 

a 

N-2 

a 



, a =exp (-j 


is the N X N DFT matrix with 'a' as the Nth root of unity. 


X = X 


( 2 . 8 ) 


is the inverse DFT of X. 


DFT of the correlation function p is: 

■^-uv 


1 

N 


N-1 


^uv^^^ = N exp(-j 

s=0 


. 2 t ks 


N 


) 


1 

N 


N-1 N-1 

S ( S ) exp(-j 

s=0 i=0 ^ ^ ® 


, 22iks 

N 


) 


Substitute i-s -t, then^ 


P,„(k) = i V* exp(-j 

i=0 ^ t=0 ^ ^ 


uv 


N 


1 

N 


4 2 7t ik v ^7,^ * , . 2jnct^ 

^Sq 635) (-J e3{p (+j ) 


N' IT, . VI 
k k 


(2.9) 


where U and ¥ are the DFTs of u and v^ respectively. 
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Hence, P = N (U . V*) 

~uv — — 

where U = ^ F u and Y* = ^ F* . 
Thus, = N F* (U . V*) 

This can be stated as a Lemma. 


( 2 . 10 } 


(2;ii) 


Lemma 2. 1 ; 

The cyclic correlation between two complex valued N- length 
sequences u and y is N times the inverse DFT of termwise product 
of U and V*. 

Corollary 2.1 : 

If u and y are real- valued, then, 

Yuv (V* = I F* y) (2.12) 

p^^ = N F* (U.V*) (2.13) 

Corollary 2.2: 

If u = V, p = p is the auto- correlation function and 

thus , 

= N ( U.U*) = N lyl ^ (2.14) 

2.1.1 Conditions on DFTs of Sequences for Ideal Correlations 
Ideal auto-correlation, from (2.4) is: 

Ea = (H O ....* 0 )'^ 

For u to possess it, it is necessary that: 



2 


IJ 


P 

— u 




F p 
= 


or |yl ^ F p 


= ~ F . (N 0 . . . 0)"^ 
= N (11..... 1)^ 


SO, |y| =(-i)^ (11 1)^ 


Lemma 2. 2 : 

A necessary and sufficient condition for a complex- valued 
N-length sequence u to have ideal cyclic auto- correlation iss 

1 ^ 

( jj arbitrary, 

i.e. each spectral component of u should have uniform 

"1 JsC 

magnitude ( ^ )'^. 

Similarly, ideal cross-correlation, from (2.5), is: 

Euv ° 

so P, = N(U.V*) 

= i F p 
N = s-uv 

= (0 0 ..... 0 )'^ 

or H‘V* “ (0 0 ..... 0)^ 
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Lemma 2.3 : 

Two N-length complex- valued sequences u and v possess 
ideal cyclic cross- correlation iff for each i, either one or 
both of and are zero. 

These results indicate the transform approach to correla- 
tion analysis. Conditions for ideal correlations have been 
specified. 

Thus/ ideal cyclic correlation sequences should possess 
the appropriate DPTs. Ideally, sequences with both ideal auto- and 
cross-correlations are desired. No known sequences' possess ; 

both these assets. A class of poly-phase sequences called the i 
Prank-Zadof f-Chu (PZC) sequences and Nth roots of unity 
sequences are known to possess ideal auto- and cross-correla- 
tions, respectively. These are discussed in the next section. i 


2.1.2 PZC and Nth Roots of Unity - Sequences 

PZC sequences for length N are defined as follows 

[chu, 197^: 


(i) If N is even and gcd (M,N) = 1, 

2 




then u. 


k 


exp (j ) 


(ii) If N is odd and gcd (M,N) = i, 
then = 63^ ( j , , 

Example 2.1 ; Pop instance# consider the i-leng-th sequence (M=l) ; 
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= 1, = e3<p(j 45®), =-l, U2=e3<p(j 45®). 

The DFT, U, of this sequence is: 

U = ih exp (j 45°) % % exp ( j 135°) h)'^ 

Also p^= ( 4 0 0 0 )^ 


In general, it can be shown |^Sarwate, 1970 for the FZC 
sequences that: 


p^ = (N 0 0 ... 0)^ and 

Uf = (^)^ e 3 ^(ja^) as required by Lemma 2.2. 


Similarly, another class of poly-phase N-length sequences 
formed using the Nth roots of unity can be seen to have ideal 
cross- correlation. There are a set of N such sequences for 
every N, Each N-length sequence is a row of the N x N DFT 
matrix. So, the tth element of the ith sequence. 


'x^ - exp (- j - "^^" ) . 


Consider two sequences u and v: 


u = 2 , the ith sequence and 


m 


V s= X , the mth sequence, 

"k 


/ . 2 w±k^ / . 2 mk 

So, u^= e5q>(-j — — ■), v^, = e3q)(-ji 


1 

U ;= ~ F 
“ N i - 


N 


N 


V 


* 1 T-i* 

-■^1 2 ‘ 


Thus, (U , *■ ( 0 


0 i-^1 

i 

i4-l 

N-l 

( 0 • • • • 0 

1 

0 * * 

• 0 ) 

0 iti- 1 

m 

m-f 1 

N-l 

( 0 * . • » 0 

1 

0 * 

» * 0 ) 

( 0 0 

0 

0 - 

••0 ) 


) 

T 

T 

T 
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This satisfies the condition 'of Lemma 2.3. Hence, these 
sequences are absolutely uncorrelated. 

Unfortunately, these poly-phase sequences are of 'little 
applicational value as they are hard to generate. So, easily 
generated binary sequences with sub-optimal coi~relation 
properties are to be accepted. LFSR sequences are the most 
widely used such sequences. In fact , their mathematical and 
operational elegance leaves little else to be desired as far 
as cyclic correlation is concerned. 


2.2 MAXIMAL- LENGTH LFSR SEQUENCES 

These sequences, also called m-sequences or pseudo-random 
binary sequences (PRBSs) , are ■(^•'1^ (+1, -i) sequences genera- 

ted using linear- feedback shift-register circuits. They are 
known to possess asymptotically ideal (near-ideal) correlation. 
Asymptotically-ideal auto- correlation is defined as follows; 


Let ^p, 


|^maxjp^(s)j : l^s^N-1^, 


max 

then u has asymptotically ideal auto-correlation if; 
a 


lim ?max 

N -»o0 p„ (o) 

u. 


0 


(2. 15) 


Before going on to their asyni> to ti call y-ideal correlation, we 
give some preliminaries about m*- sequences. 



2.2.1 Generation of m- Sequences 


The general configuration of an r~stage 2 -1 -length 
sequence generator LFSR circuit is as shown in Fig. 2. 1. 



J / 


This circuit generates the sequence (%}' 


n — 0/ 1/ 2/ 


with periodicity (2 -!)• Output u^ for an r-stage LFSR is 
given by the following recurrence relation as is evident from 
the configuration is the exclusive OR summation) : 


% ® "=i %-i 

1=1 


0 ; u^r/ c^G^O/ 1^ (2.16) 


Define operator D as: Du^ = ^i+l ' 

Then (2.16) becomes: 

(D^ ® ^ D^"'^)u^ = 0/ n^O 


i=l 


(2.17) 


The solution of this I^Key, is obtained using the extended 

field GF (2^) and the initial loading (u^u^. . . . 

From (2.17i^the char aeteris tic equation of the generator 
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X 


r 

S 

i=l 


X 


r-i 


0 


(2. 18) 


This is also the connection polynomial which specifies the 
configuration of the LFSR circuit, Alongwith the initial 
loading, it specifies the generated sequence completely. Initial 
state has to be a non-zero state to obtain a non-trivial solu- 
tion. LFSR being a finite- state m/c, the sequence of its states 
has to be periodic. The maximum period of an r-stage LFSR can 
be 2-1 as the all- zero state is excluded. Hence, the name 
maximal- length sequences for such outputs. This shows that 
not all connecting polynomials generate m-seq\iences. In fact, 
a necessary and sufficient condition for m-sequence generation 
is that the connecting polynomial should be a primitive irre- 
ducible polynomial of GF (2^) .Qsolomb, 1967]. To illustrate 
this, we consider two exanples - one using a primitive poly- 
nomial and the other a non-primitive irreducible one. 


Example 2.2 : Consider the non-primitive irreducible polynomial: 

x^ ® x^ ® x^ ® X ® 1 = 0 

Initial loading is (1 0 0 O) . 

The LFSR configuration is shown in Pig. 2.2(a) and the 
sequence generated is: 

u u 
o 

0 0 0 1 1 0 0 0 1 1 . . . 

Period of the sequence is 5 and not 2^-1 - 15, 
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F!G 2*2(a) 



Fig 2*2,{jb) 







Example 2*3; 

The configuration of Fig, 2.2(b) is for the primitive 
, 4< 

polynomial of GF(2 )s 

® X ® 1 a= 0. 

Initial loading is ( 1 0 0 0 ) . 

The generated sequence is 

%^r*- ^14 

000 100 110 101111000 1 .. 

The period obviously is 15. 

This brings us to the problem of finding the number of 
distinct (cyclic shifts of a sequeice belong to the same equi- 
valence class) m-sequences of length 2-1. This is a well knowi'^ 
result. 


L^tuma 2.4 fWlomb# 1967|4 


The number of distinct mi-sequences of period 2 -1 iss 

g((2^-l) 

r ' 


where ^(n) is the Ealer^totient function defined ass 



if na*l 

p^i-i If n^ 1 


■and by the' uni^que f actbrls'itioti theorem,', eveiry integer n"^- 1 
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n = n 

i=l ^ 

JZI'(2^-l) is the number of positive integers less than 

JT 

2-1 and prime to it. 

0’(2^-l) = 2^-2 if 2^-1 is a prime. 

(called Mersenne prime) . 

Hence, for prime periods 2^-1, the number of distinct m-sequences 
is: 


r 

With these preliminaries of m-sequence generation, we go on to 
enumerate some of their important properties. 

2.2.2 Properties of m- Sequences 

Some of their important properties exhibiting their 
mathematical elegance and potential for use are: 

PROPERTY l : SHIFT PROPERTY 

Cycle shift T®u of sequence u is another m-sequence. 

PROPERTY 2 ; SHIFT AND ADD PROPERTY 

Mod 2 sum u ® T®u of a sequence and its cyclic shift is 
another "phase” of u. 

PROPERTY 3 t WINDOW PROPERTY 

An r-length window moved through a (2^-1) -length 
m-sequence ' shows ' every non-zerp r- tuple exactly once. 



Alternatively, each non-zero state occurs exactly once in a 
period of the sequences 


e.g. u = 000 100 1 10 10 1111000 1... 
of exarcpie 2.3 exhibits the following sequence of states; 


0 

0 

0 

1 

1 0 

1 

1 

0 

! ° 

.1 

1 

1 

0 

0 

1 

0 

4 

1 

1 

0 

1 

4 

1 

1 

1 

1 

0 

1 

0 

0 

1 

0 

1 

0 

1 

1 

1 

0 

1 

0 

0 

1 

0 

1 

0 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

0 

1 

1 

1 

0 

0 

0 


PROPERTY 4: AUTO- CORRELATION PROPERTY 


Consider the 


{,. .) 




m- sequence u formed out of the 


u using: u. 


1-2 u^. 


Then, the auto- correlation of u is: 


P c^s) 


I N, s = 0 

^“1, S — 1/ 2/ a « . , (N— l) a 


r, a 

SO, 


1 and 


•^itax _ lim 1 
N ^ CO P^(o J ^ N 


Hence, from (2.15), these sequences have near-ideal auto- 
correlation. This can also be seen in terms of their DPT. 

For the given p^(s). 


P, tk) = i i p, .( 1 ) exp(-j iril) 


(a) P^(}t»0) 


i*0 

1 N-1 

# S P 
N . " ^u 


N (N-(N-l)) = I 



2 it ikv » 
1 ) 
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(b) p^(Mo) 



N-1 

2 exp{-j 
i=l 


So, (Mo) 

Now, P^ (k) 


■ N-1 

S exp (- j 
i=l 

I (N -(-!)) =1 

NjU^l^ 


2 ^ik 

n" 


) 


(N + 1) 


Thus, juj^ = 


-4 , k=:0 mod N 

N'^ 

^ (N + 1) ,, k/^0 mod N 
N^ 


- 1 for k;»^0 . 


(2. 19) 


Also, all sequences having (2.19) as the DFT possess the near- 
ideal correlation. This leads to the Lemma; 


Lemma 2.5 ; 

The periodic auto-correlation of a sequence u is of the 
m- sequence kind iff; 



1 

H exp ( j a^) / k=0 mod N 

(a^ arbitrary) 

^ (N+l)^ exp(ja^), kj^O mod N 


Example 2.4 ; Let u = (-1 1 
TT /I J180“ 2f2 j249® 2f2 

vJ '^1 # *"**"7^'*' 


1 -1 1 -1 -1) , N =7, 
j249® 2f2 Jill* 2 I 2 

zM Ml 

76,7 


249'* 

# 

111 “) 


This transform approach to m-sequences is not popular in the 
literature.' As is ''evident from Lemma 2,5, this gives a generalised 
approach to im- sequence- type correlation sequences. Non-binary 
sequences have also been incorporated in the general results of 
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designing generalised cyclic sequences with desired ideal and 
near ideal correlations. 

This auto- correlation property is useful for many applications 
as mentioned in the first chapter. For instance/ in Radar 
Ranging/ a long m- sequence (depending upon the maocimiom range 
to be estimated) is transmitted. The received sequence is a 
delayed version of the original/ delay being directly propor- 
tional to the range of the target. This is correlated with a 
replica of the original and time-shifts counted till a peak is 
obtained. This gives an estimate of the target range. Another 
potential use is in the identification of linear cyclic-shift 
invariant (LCSIV) systems. 


2.3 LCSIV SYSTEMS 



Let Z be an N-point LCSIV system. Then, 


1 

w. = a 
^ J=0 

^i-J ^j , 

N-1 ( 

= ' S 

^i-j J 

j=0 

Alternatively# 

=i 'Z . X 'r 


/ (i-j) mod N. 


( 2 . 20 ) 


■'ij'^N xN ' 


Zij = <f (i-j) mod N. 


Rows of Z are Gyclic shifts of each other. 
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For instance/ for N = 4 


^o 


— 

(SI 

0 

^3 

^2 

^1' 



^ 1 : 

: 


^0 


^2 : 



^2 


^2 



^3; 

' 

^ 2 j 

■ X3 

j 

"3 

^2 

Zi 

^0; 


W 3 :i 


2.3.1 Elementary Properties of LCSIV Systems 

PROPERTY 1 : DFT row and column vectors form an eigensignal set 
for LCSIV systems. 

PROOF: Consider V = F Z F~^ / where F is the N x N DFT matrix, 

F F* 

SO/ Y 
Let G 
so, 


N I 

I I i l’ 


I I 

N-1 
S 2 


m=0 


f* 
im mk 


^nk 


1 

“ N 

N-1 

S 

i=0 

^ni 

%k 


II 

N-1 

S 

i=0 

fni 

N-l 

2 

ms=0 

z. 

xm 

II 

N-1 

S 

i=0 

N-l 

£ 

m=0 

^ni 

^mk 

II 

N-1 

s 

i=0 

N-l 

2 

msO 

^ni 

* 

^mk 

1 

N-l 

N-l 




Z-. . 


2 . 


(f* f ) / (f f =1 and 
nm nm^ ' nm nm 


2 


im 


2# / wn ^ 


« (i-m) ^i~m (k-n) 

iasO m**u 
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Yr 




N-1 



= 

1 

N 

S 

m=0 

m. 

(k-n) 



Let 

i-m 

= s. 


1 

N 

N-1 




S 

m=0 

m. 

(k-n) 



N-1 



rr 


S 

m=0 

f* 

m. 

(k-n) 

Z 

is 

the 

DFT 

of (! 



f 

2n 

, n=k 

^nk ~ 

1 0 

, n;^k 


N-1 
2 f 
i=0 


n, (i-m) ^i-m 


N-1 

E f z 
„ ns s 
s=0 


/ (Z, 




N-1 

n " ^ns ^ 

s=u 


-N-l' 


Therefore, ^ Y is a diagonal matrix with diagonal entries as 

the DFT values of the LCSIV systom vector z. The similarity 

transformation F Z F uses F to diagonalize Z. Hence, the 

rows and columns of F are eigenvectors of Z. 

~ Q.E.D. 


PROPERTY 2 : w = N F* (Z . X) 

i»e, the output w is N times the inverse DFT of the termwise 
product of DFTs of the system vector z and input x* 


This follows easily from 1. 

property 3: If O^ is the cyclic-shift operator i.e., 

S Hi “ * (i+h) mod N 

which is ,a cyclic- shift of sequence u, the system matrix Z 
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comrnutes with 0 

I 0 X = 0 I X 

This is the cyclic-shift invariant property and follows readily 
from (2. 20) . 
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= s £ x._j x..^ 


2 z. 

j J 


X is the perturbation signal to be chosen to have ideal or 
near-ideal auto- correlation. Two cVioices .are of interest: 

(i) Let X be an N-length FZC- sequence 
then, p^(j-k) = N 6^^ 


So, e, = N z. 


k ' 




N 


Each is proportional to z-^. 

(ii) Choose x to be an N-length m-- sequence (N = 2^-l) 


then, ^ ^ 2j 


j/k 


In other words > 


— 






N 

-X . . . . -1 : 


^o 


®o 

-1 

N .... -1 




®X 

» 






• 





# 

: -1 

' X * ♦ ' • * , H' 




f 

f ' ' • 

i ' ■ 

1 ®N-X, 

. J 


Solution of this 

3.:, d 


Zq 






'' * 

^N41; 



UMh . : 

i . 't 


2/ (N^l) 1/ CU+l) . . , X/ (N+l) 

1/ (N+ 1) 2/ (N+l> , . . 1/ (N+1) 

' ♦ • 

1/ (N+ X) 1/ (N+1) , . . 2/ (N+1) 


o 


'N-X 
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i.e. - nIt + N+1 ®j 

So, z = T e 

The interrelation of sequences with good cyclic correlation 
properties, LCSIV systems and DFT has been presented here as 
a prelude to a similar study of LDSIV systems and sequences with 
good dyadic correlation properties. Before we go on to that in 
the following chapters, we present briefly some sequences with 
good aperiodic correlation properties. 


2.4 SEQUENCES WITH GOOD APERIODIC CORRELATION PROPERTIES 

Recent advances in spread- spectrum code-division multiple 
access (SS-CDMA) communications have focussed attention on the 
need for large sets of-|^+l, -1^ -valued sequences having good 
periodic and aperiodic correlation properties to be used as 
signatures for various users. A detailed analysis of the 
requirements for such sequences and inportant correlation para- 
meters of LFSR and relate<2 sequences are presented in |^arwate 
and Pursley, 1980 and |0tirsley and Roefs, 197^ . Some salient 
points are presented here. 

-?!periodic auto- and' cross-correlation of u = (u^. ..Uj^_^) 
and V = (v^,.,v„ ^) are defined as: 


C^(s)^ 



N-1 


: 0 , fsf^N . 
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Ideal and near-ideal correlations are defined similar to the 
periodic case, 

-1^ Barker sequences are knowi to have good auto- 
correlation. But they are not known to exist beyond length 13. 
If u is a Barker sequence then; 

f N , s =0 

(s) = < 0 , s odd 

1 , s even 


Example 2,5 ; The length- 13 Barker sequence, u and are : 


u 1 1 1 1 1-1-1 1 1-1 1-11 

C^(s) 13 0 10 10 10 10 10 1 

Generalized Barker sequences have also been studied but 


are of limited practical value. 


Another class of aperiodic sequences is Huffman sequences 
If U is such an M-length sequence then: 


■ C 


, s«0 

/S — 1/ (N— 2 ) 

, S=(N-l) 
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Example 2.6 ; A length-5 integer Huffman sequence u is: 

u 1 2 2 -2 1 

14 0 0 0 1 

f 

But these sequences are non-binary ones and hence not used 
frequently. 

The choice again falls on LFSR and related sequences like 
Gold and Kasami seqa.ences. These sequences have, the near ideal, 
cross- correlation property: 

lim max n ^ 

^ ° • 

Large sets of these sequences are known as N becomes larger and 
larger. Hence they are suitable for SS-CDMA systems. Tables 
of correlation parameters for LPSR/ Gold and Kasami sequences 
upto length 255(2^-1) are given in |^Pursley and Roefs, 1979^- 

This brings us to an end of the background material on 
sequences with good periodic and aperiodic correlation proper- 
ties. Now we go on to the domain of dyadic correlation and 
LDSIV systems and analyse a class of Boolean functions for their 
dyadic correlation properties. 



CHAPTER 3 


DYADIC CORRELATION AND LDSIV SYSTEMS 


This chapter deals with dyadic correlation/ LDSIV systems 
and their relation to Walsh-Hadamard transform (WHT) . After 
presenting dyadic correlation and its relation to LDSIV systems/ 
in Section 3.1/ we review Walsh-Hadamard transform in Section 
3.2. Elementary properties of dyadic correlation and LDSIV 
systems in terms of WHT are ej^ounded in Section 3.3. Section 
3.4 gives a scheme for identification of LDSIV systems using 
dyadic correlation and sequences with ideal dyadic auto-correla- 
tion property. A typical hardware realisation of a 64-bit 
dyadic correlator constitutes Section 3.5, Finally/ Section 3.6 
presents briefly the WHT-based classification of Boolean func- 
tions, which is to prove useful in dyadic correlation analysis 
of NRQFs in the next ch^ter. 


3.1 


DYADIC CORRELATION 

Dyadic correlatibn/ between two N-point (N=2^) 


sequences ^ 


V 




) 


*'1 


is defined as: 
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1 


E u . V . « 
1=0 ^ 


(3.1) 


where i = (i^ ij^_^.... i^) , s = (s^ Sj^_^.... s^) are the 
binary representations of i and S/ 

les = (i^es^ ijSsi), 


(i„®s„ i^ 1 ®s„ . 
n n n— 1 n- 1 


i^®s^) are the bitwise 


mod 2 difference and sum of i and s; 6 and ® are same mod 2, 


In terms of the dyadic matrix, V, 




Example 3,1 ; Let u = (u^u^u^u^^ Z = *'^o’^l'^2^3^ 

4-length sequences. 


■^o ^2 '^3 

'"o ^3 ^"2 

v 

2 3 o 1 

V „ V V 

3 2 1 o 


Similar to the cyclic case, ideal dyadic auto- correlation for 
an N-lmgth sequence is defined ass 


h, (s) « N 6 , s=0, l, . . (N-1) 

u ' 'SO 


(3. 2) 


Ideal croi^S‘-<'COirrelati' on," is,* 


h^^ Cs) *** 0 / ssO , 1, • • t f (N— l) 


(3.3) 
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Sequences with such properties are to be discussed later. 

Dyadic correlation# which is shown to be the same as 
dyadic convolution, relates the output to the input of linear 
dyadic— shift invariant (LDSIV) systems. These systems are 
characterised by dyadic matrices. If Y is an N x N (N=2^) 
LDSIV system, then. 


u 



> V 


(CONVOLUTION) 


(CORRELATION) 


N-1 

V. =: £ V. ^ . U . 

^ j=0 

N-1 

= S Y.q, u. 

j=0 J 

N-1 

Alternatively, v,= Y u, where Y = (y. .) 


ij N X N ' ^ij 




Hence# properties of LDSIV syst^as ar^ nothing but those 
of the operator Y which also defines dyadic correlation. 

Various aspects of dyadic correlation and LDSIV systems 
can be studied using Wal sh-Hadamard transform which is reviewed 
here briefly. 


3.2 WALSH-HADAMARD TRANSFORM 

The Walsh-^Hadamerd Traasfo (WHT) ^Ahmed and Rao, 1970 

Jr, 'Itl' ' 

of aa N a 2 ;^point signal represents its seguency-content or the 




Walsh functions' content among the first N i.e. 0 to N-1 Walsh 
functions. Each row of the transfoimn matrix is an N-point 
Walsh function. This is similar to the frequency- spectrum, 
represented by the Fourier transform, of a signal. WHT can be 
defined with respect to the following three different orderings 
of the Walsh- functions in the transform matrix: 

( l) Sequency- Ordered Transform 

Rows of this transform matrix, W, are Walsh- functions in 
ascending order of their number of zero-crossings i.e. the 
sequency. The ith row W^^ is the Walsh function with i zero- 
crossings. For an N-point (N=2^) signal, X' this transform, 

X^, is: 

2^ = ^ ~ ' where W is the N x N sequency- 

ordered transform matrix defined below: 




(- 1 ) 


n-1 
k=0 


^^n-k ® Vk-l^J’k 


(3.4) 


where i and j have the binary representation: 


i= 3 (i n.....i,r},iL 

' n-1 n-2 *••*• -^1 o'' 


n -1 , 

S i,.2^ 

k.O 


^ ^n-2 


n-1 V 
2 V2 , 


ik^Jk^tO, ij 


and 0 . 


A few ' trivi.aA ;|5'ropeftie's. of W are: 


::^i3 


(symmetry) 


(3.5) 
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(ii) 

"'^i j * 

^ik 

- ^i,j9^k 


(CLOSURE) 

(3-6) 


N-1 







(iii) 

2 

k=0 

^ik 

w = N 6 . . 
Jb ij 


(ORTHOGONALITY) 

(3.7) 


or 

1) 

« 

SH 

I 





or 


1 

“ N 

W 



(3.8) 


So, 

X = 

= M 2 w 

- 




Exantple 3. 


The 4 ; 

X 4 matrix 

W with its sequency 

ordering is 







Sequency 




1 

1 1 

1 

’ Ho 

0 



W = 

1 

1 -1 

-1 

Hi 

1 




1 

-1 , -1 

1 

H2 

2 



i 

1 

1 

-1 1 

~1 

Ha 

3 



(2) Pal ey- Ordered Transform 

The transform matrix P results in the transform: 

”^P ^ ^ ^ 


where 



(- 1 ) 


n-1 

S 

k=0 




(3.9) 


A oomp arisen of (3,4) , and (3,9) shows that if 
t - (t t ^... t.t_) such that t = ® ' s=0/ l/..(n-l) 

then P^,. 

But t is the Gray code of i denoted by i^. 


'.'Thus,/ Wj 


-i 


G 
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Conversely-/ i = t © i for -the given -t; s=0/ 1/ . . . (n-l). 
s s s+1 

This is the reversed Gray code/ ^t^ , of t. 

Thus/ P+. = W . 

—■c r , 

G 

Example 3.3 ; For N =4, 








i 

r . 

^G 



1 

1 

1 

1 

£o 

0=00 

00=0 

Ho 

p = 

1 

1 

-1 

-1 

£l 

1=01 

01=1 

,Hi 


1 

-1 

1 

-1 1 

^2 

2=10 

11=3 

Hs 


1 

-1 

-1 

1 

^3 

3=11 

10 = 2 

H2 


P also satisfies the proper-ties (3.5)/ (3.6), (3.7) and (3.8), 


(3) Hadamard-Ordered Transform 

This transform matrix, H, has proper-ties similar to W 

and P/ 


X, 




^ H X / where, h^^^j = (- 

n-l 


(- 1 ) 


k=0 


<^/ is the inner product mod 2 of i and J. 

(3.9) and (3.10)* it is seen that if 


ia == Xi l..*.i .)/ the Mt^reversed i, then, 

'' 'R ' ' o ' 1 








(3, 10) 
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Example 3,4 ; For N = 4 , 



Rows and columns of each of the transform matrices are 
linear Boolean functions (LBFs) to be defined at a later point 
in this chapter. In addition^ rows of H are in terms of the 
inner product of indices which is a natural representation for 
LBFs. This property of Hadamard-ordered transform has been 
used throughout this work. 

Having dealt with the preliminaries of various transforms/ 
we pass on to an analysis of dyadic correlation and LDSIV 
systems using the Hadamard transform. 

3.3 ELEMENTARY PROPERTIES OF DYADIC CORRELATION AND LDSIV 
SYSTEMS 

PROPERTY 1 8 EIGSNSIGNALS 

N-poinfe Walsh functions are the eigensignals of LDSIV 
syst^s# i.e. the rows and columns of Walsh matrices are eigen- 
vectors of dyadic matriees. 

PROOF ; Consider I J , where H is the N x N Hadamard 

matrix and Y is the system 


matrix. 
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G 


I H Y H (as i 


aik 


II 

N-l 


N-l 





S 

i=0 

h . 
ni 

E 

m=0 

Y- 

-*^xm 

^mk 



1 

N-l 

N-l 






X 

- N 

S 

i=0 

E 

m=0 

^im 

h , 
nx 

^mk 



■1 

N-l 

N-l 




9 


1 

“ N 

E 

i=0 

E 

m=0 

^im 

h , 

nx 

\k 

Km 

= 1) 


N-l 

N-l 






X 

““ N 

E 

i=0 

S 

m=0 

"'^n^ i®m 

^i®m 

^m,n®k ' ^^im ~ 

^i®m 


m=0 


h 


N-l 

S 


m,n®}c “n,i®m ^i®m 


CLOSURE property 
of H ) . 


Let i®m 


Sf then, 
. N-l 


N-t1 


< 3 . 


nk 


N ^m,n®k \s 


N-l ' 

^m,n®k ' 

m=0 


1 

where is the nth Hadamard coefficient ofyasY = |jHy, 
^ ss • 'Yn-I^ system sanple vector. 




N-l 

Ihus, 

%k 

m^n 

m=o 


N-l 

h a,, =s n 6 , , because sum of any row or 

" in./ .iriiwiK nic 

Now 

S 

m-0 


ooixinm entries of m H matrix is zero except the zeroth row and 
column. 



So, g. = N Y 6 T .. 
^nk n nk 




Yq 0 .... 0 


0 Y^ 0 . 


0 


0 0 0 ... Y, 


N-1 


G is a diagonal matrix with the diagonal entries as the WHT- 

values of the system vector y. We have been able to diagonalize 

1 


Y using., the similarity 'transformation H Y H 
and columns of H are eigenvectors of Y. 


Hence, the rov.rs 


Q.E.D. 


Example 3.5 ; Let Y be a 4 x 4 LDSIV system matrix .. Choosing 
the third row of H, , we get 

r - 
1 


Y H 


T 


Yo Yi 72 Y 3 
Y^ Yq Y3 Y2 
^2 ^3 ^o 


-I 

-1 

1 


= (Y0-Y1-Y2+Y3) 


1 

-1 

-1 

1 


So, is one of the eigenvectors with the eigenvalue 
' Also, as a verification: 


H Y H“'^ 


0 

0 


0 

o 


yo-Yl+y 2 ”y 3 ° 

o yo+yry2-y3 o 

o 0 Y^Yi-y2+Y3 


O 


0 
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The ith diagonal entry is the ith Hadamard coefficient 

of 2 = (^0^1^223^ 

PROPERTY 2 ; DYADIC- SHIFT INVARIANT PROPERTY 

)c . k 

If R is the dyadic shift operator such that R = 

— i©k the kth dyadic shift of u, then the dyadic 

correlation operator, Y, commutes with R , 

Y r’^ u = R^ Y u . 


^iej ’ 


Q.E.D, 

PROPERTY 3 : TRANSFORM : REPRESENTATION OP DYADIC CORRELATION 
Ttie dyadie correlation, v, between u and j is N times 
t^o' inverse: WHT , of WHTs, U and Y, of u and y . 

i.e. V = N H (y.u) , U = i H u , 

i = i S ^ , 


PROOF: 


V = Y u 


Let = Y R u 


N-1 


Thus, ' 


Let j®k = t, so, j =t®k . 


Thus, 


N-1 

2 y. 


t=o 


i®t®k 


Ut , (y, 


ij 


(r"^ y u)^ 


Hence, Y R^ u #= R^ Y u 





This property is evident from property 1, As Walsh 
functions form the eigensignal set for Y , so the input-output 
of an LDSIV system are related through the WHT similar to LTIV 
systams where sinusoids are the eigensignals and, accordingly, 
Fourier transform relates the input-output. 

Now ’((^e go on to detponstrate the use of ideal dyadic auto- 


correlation sequences in the identification of LDSIV systems. 
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3.4 IDENTIFICATION OF LDSIV SYSTEMS 

A block schematic for this purpose is shown in Fig. 3.,1, 
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b (s) = N 6 ^ / then 
u so 

= N 

Thus, y^= i 


In terms of the transformation T, 



Hence/ all components of y are coirpletely Icnown. 

The identification requires a dyadic correlator. A 
typical circuit realization of a 64-bit dyadic correlator is 
presented next. 


3,5 THE DYADIC CORREDATOR 

Dyadic correlation between two N = 2^-polnt functions f 
and 2 can be irrplemented in two ways: 


(i) ComputationallY using the fast Walsh-Hadamard Transfo3rm(PWHT) | 


Ahmed and Rao 


, 1975^ 


, ' ' ''.I ' r' ■■■ '■ ^ 

s 5 ='N s 2 


E/ g & H can be inplemented using the FV3HT very 


efficiently. 
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(ii) Using digital hardware to implement the hard-wired 
cori/elator. This uses the definition: 

b. (s) = "s'" (.i,f(x)®g(x®s) 

^9 x=0 

ra (Number of zeros in f (x)<Sg (x®s) ) - 
(number of ones in f(x)®g(x®s)) 

A t;ypical circuit for 2 -point functions using ‘16-1 MUXSS, shift- 
registers and counters is shGiy’n in Fig. 3., 2. 

The correlator consist3 of a PRE-SET REGISTER which 
contains s = (s^s^. . . . Sg) / the dyadic shift for which the 
correlation is desired. These bits are X-ORBD with the output 
of a mod-64 counter counting from Q thru 63 to gi'^e the six 1 

control bits c^ to c^ for the MUXES. c^ and c^ control the | 

1-4 DEMUX whose outputs are the strobe inputs of th« 4 16-1 i: 

MUXES. STj,=0 means the MUX is disabled and ST^=:1 enables the | 

MUX. ST^ along with Cg to Cg connect the output of the suitable i 
MUX to g (x®s) as X goes from 0 thru 63. Each g (x®s) is XDR~ed | 
successively with f^ to This output governs the count of j 

a 7-bit up-^down counter pre-set at 63. A 'O' counts up and a 
' 1* counts down. Hence, the output of the counter after 64 clock 
pulses IS the difference be twe^ number of zeros and ones of | 

' ' ' ' ' ,.v ' ' ' ' ' ' \ 

(g(x®s) + f(x))plt^ the preset number 63. . Subtracting 63 out | 

of it gives tile dj^'idiC b^g(s). f 

^ correlator for any 2^-point functions can be i 

realizsedf I 
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NOW/ WHT-based classification of binary Boolean functions 
is briefly dealt with as it will be of use in dyadic correlation 
analysis later. 


3.6 NOTE ON THE CLASSIFICATION OF BOOLEAN FUNCTIONS 

Boolean functions can be divided into various equivalence 
classes in terms of their Walsh-Hadamard Transforms |^Karpovsky/ 
197 sj. As mentioned earlier/ each row of H, , is the ith 
linear binary Boolean function in the (i. -i] - form. The cth 
linear Boolean function of n-variables in |0/ -form is 
defined as: 


® (c) / \ 

f (x) = c. 

i=0 


(3. 11) 


where c = (c„ c„ ^ 
n- 1 n- 2 


s’- 


X 


= ^-2 ^o^ ^ ^i' x=0/ l/.../(2^-l) 


AlsO/ (x) =r<^ C/ X 


Its 


{l, -l} 


form isj 
(c) 


How, ^ 


Inch hij = (-1) 




n 


SO/ as j yaries from 0 to 2 -1/ 


(x) , the ith LBP. 



For any N( 2^) -point Boolean function g, 
2 ~ N = — 


G. 


1 


1 

N 


N-1 

S 

j=0 



i — 0 / 1 / 


, (N- 1) . 


Thus, each is the normalised correlation betaveen the ith 
LBF in -ij- -form and g. It also gives the distance between 

g and the ith LBF. These G's can be used to classify Boolean 
functions . 


Example .3.6 ; Consider the classification of LBFs themselves 
LBFs will be classified both in their -[o, l]- and fonns. 

The ith LBF can be written as: 


®f = 1 (H^ -H^) 

(i) 1 IT (i) 

SO, F = - H . t 


(3. 12) 


®T- (i) 

^ Jt 




I • I ( <Sj ' iO -<Sj ' Si>^ 


-i. (N ^ . -N ^ ) 

2N ^ jo Ji 


4 ( 6 . -6.,) 

2 ^ jo ji 


Kenee> ' ''WHI of the itfe-I^BF is: 


® (i) 

j 


1/2 , i =0 

0 , , j/0 , i 7 j ? 1# 2, . 


(N-1) 


(3.13) 
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This can be e5<pressed as; 


Lemma 3.1: 


A binary function f in the 


1/2 , 3=0 




l\ -form is an LBP iff: 
J 


-1/2 , for some j =k 

0 / j /O, k , 


where P 


Then it is the kth LBP. 


.e 3.7; Consider the 3- variable 7th- LBP, 
®f (x) = x^® Xj^® X 2 / k = (111) =7. 


000 001 010 011 100 101 110 111 


®£( 7 ) 0 1 


0 10 0 


®p(7) 


1/2 0 0 0 0 0 0 - 1/2 


Classification of LBPs in the 
evident. The ith LBP is; 


1^+1, -1^ 


form is even more 


tA(l) 

J 

so, 




N-1 

t=0 


Thus, P, 


<i) 


N * ^ ^ji 

^ ji 


1 / 3 i 


0 , otherwise 


(3. 14) 
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Lemma 3.2: 


A binary function f in the ^+1/ 


-1^-form is the kth LBF 


iff: 


j 1 , j =k 

1^0 , otherwise. 


Example 3.8 ; Another class of ^+1, -l^V functions called Bent 
functions are classified with respect to their WHT, 


{■^ 1 ' - 


is a Bent function of n variables iff; 




1 2"-i. 


f(x) = X^X 2 ®X 3 X 4 ®X^X 3 ®X^X^®X 2 X 3 ®X 2 X 4 is one of the bent 
functions of four variables. 

This brings to an end the exposition of WHT techniques 
and their application to correlation analysis. Now we go on 
to use these for a detailed dyadic correlation analysis of 
non- r^etitive quadratic forms. 



CHAPTER 4 


NON- REPETITIVE QUADRATIC FORMS 

This chapter is devoted to an analysis of dyadic correla- 
tion and related properties of non- repetitive quadratic forms 
(NRQFs) . in Section 4» 1 we define NRQFs and give some prelimi- 
nary details about them - viz. principal NRQFs (p-NRQFs) . 
special NRQF (s-NRQF) . elementary properties, generation of 
NRQFs using counters and minimal LFSR circuits. Walsh Hadamard 
transforms of NRQFs are computed in Section 4.2. In Section 4.^ 
these are used to derive dyadic auto- correlation of NRQFs. 

Their cross-correlation properties are abounded in Section 4.4 
using permutation cycle structure relationship^ between t 
Section 4.5 gives the number of NRQFs related through L-length 
permutation cycles# Finally/ some results on periodic 
^eriodicS correlatibn parame^^ of p-NRQFs of 4,6 and 8 
variables are presented in Section 4.6, 


4.1 DEFINITION 

' ^ ' Non-repetitiVe <iu.a^^^ forms (NRQFs) are a special class 
of degree 2 binary Boolean functions. They are defined for an 
even number of variables# An NRQF f(x)^ of n = 2f 







expressed as a 


polynomial in the following way [Karpovshy, x976] 


f (x) 


where x 


i / j 


. 1-1 

2 x. 2 / 

i=l 


(4., 1) 


X = (x 23 ^2s-V representation or x. 

i,JS-(l, 2 2s}, a^,a.,x^,Xj6^0,l) , 


(X,) 


ai 


X. if s..= 1 
X 1 


xT if a.= 0 
i ^ 


and each term is a quadrate ^ 


Each X. occurs exactly once 
term, hence, the name non-repetltive quaaratic forms, 
the exclusive - OR adaition operation. 


.IS 


4.1.1 Principal NRQF s 

If a^ and aj are 1 for 
then the resulting NRQF is: 


all i/j€\^l' 2 ,.../n/ in ( 4 .l) 


f(x) « , 1 ^ * 

is^j 

' lr3 

such NROfs with none of the variahles conpiemanted are called 
Principal NRQPs » 

^ ^ a principal KRQF is 
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and g(x) = ®X 2 X 3 is a non 

ane talovilatsd in Tabls 4.1* 

Table 4.1 


-principal NRQF (np-NRQF) . Both 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


X3 X2 Xj) 


0 0 0 0 
0 0 0 1 
0 0 1 


0 0 


0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 


0 

1 

0 

1 

0 

1 


1 
0 
0 
1 
1 

0 0 
0 0 1 
0 10 
oil 
1 0 0 
10 1 
1 10 


f (x) 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

1 

1 

1 

0 


g (x) 

1 

1 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

1 

0 

1 

0 


IjgOTaa 4, 1 : 

Each principal form is related to the other by a permuta- 
tion of variables. 
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Hence, the principal forms are an ecjuivalence class of 
NRuFs with permutation of variables as the group operation. 
Permutations of variables can be represented using cycle 
structures or binary permutation matrices. 

Example 4. 2 t Consider the two principal NRQFs: 
f^(x) = x^X 2 ®X 2 X^®X 2 Xg and 

f 2 (x) = x^Xg®X2X2®XgX^ 

f 2 (x) is a permutation of f ^ (x) , cycle structure being: 

(l) (2 4 6) (3) (5) 


This means that x^, Xy x^ are unchanged whereas x^ is 

2 — f 4 

replaced by X 2 / Xg by x^ and X 2 l?y Xg or 
Alternatively, 

f 2 (x) = f ^ (x A) , X = (x^X 2 X 2 X^XgXg) 


where the permutation matrix is: 

r 1 0 0 0 0 0 


0 0 0 1 0 0 

0 0 1 0 0 0 

0 0 0 0 0 1 

0 0 0 0 1,0 

0 / 1 0 0 0 0 


Hence, two n variable principal NRQFs f(x) and g(x) are related 
thrcxigh the: ihvertible linear transformation: 

f(x) « g(x A), where A is an n xn permutation matrix. 
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Lemma 4.2 ; 

Each non-principal NRQF is some dyadic shifted version 
of its respective principal NRQF. 

A non-principal NRQF g (x) of n variables and its principal 
NRQF f (x) are related as: • 

where a= (a^a .....a.) and 
n- n- 1 1 

0 , if Xj^ is complemented in g(x) 

^i ~ A 

I 1 , otherwise 

In other words, g(x) = f(x®a) which is nothing but the 
ath dyadic shift of f (x) . 

In vector notation, 

g( 2 §) = f(x®a) where a is a binary n- tuple. 


Example 4.3 : Consider the p-NRQP of n = 4 variables. 


f (x) 


X 1 X 2 ®X3X4 


It generates 2' ** 16 NRQFs which are its dyadic shifts for 
a =0, 1, 15. Fifteen of .these are np-NRQPs. All of these 

in their pol,yn<x?nai as w^ array form are given in Table 4.2. 

Lemmas 4.1 and 4, 2 can be combined in: 


bemma 4.''3fe 


.'Any 'twO' 
'.'following 


NRQFs, f(x) and g(x) are related through 
of input variables: 
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Table 4,2 

A p-NRQF AND ITS DYADIC SHIFTS 

f(x) = 9x3X4 


a 5 

L-- 

(a 

•X 

3a 

f 

=1 ) I 

2 "1 1 

X 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

j f (x) 

1 

1 

1 . , 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

1 

1 

1 

0 


1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

1 

1 

0 

1 


2 

0 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

1 

1 

^1^2®^3^4 

3 

0 

0 

1 

1 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

1 

1 

5i^2®^3^4 

4 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

1 

1 

1 

0 

Q 

0 

0 

1 

^1^2®^3^4 

5 

0 

1 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

1 

1 

0 

1 

0 

0 

1 

0 

x^X2©X3X_^ 

6 

0 

1 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 

0 

0 

^1^2®^3='4 

7 

0 

1 

1 

1 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

^1^2®^3==4 

8 

1 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

x^X2®X3X^ 

9 

1 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

^1^2®^3^4 

10 

1 

0 

1 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 

0 

0 

0 

1 

0 

0 

V2®="3^4 

11 

1 

0 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

1 

0 

0 

0 

x^X2^X3x;^ 

12 

1 

1 

0 

0 

1 

1 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

^1^2®^3^4 

13 

i 

1 

0 

1 

1 

1 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

x^X2^X2S;^ 

14 

1 

1 

1 

0 

1 

0 

1 

1 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

x^X 2®^3^4 

15 

1 

1 

1 

1 . 

0 

1 

1 

1 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

^1^2®^3^4 
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g (x) = f (x A © a) , 

where A is an n x n binary invertible permutation matrix and 
a is a 1 X n binary vector. 


Example 4,4 ; The two NRQFs 

f(x) = g(x.) = x^x^®X 2 X 2 are related as: 

g (x) = f (x A © a) / where x = (x^X 2 X 2 X^) 


10 0 0 


A = 


0 

0 


0 0 1 
0 10 


and a 


0 10 0 


(ill O) 


The transformation of this Lemma just rearranges the 0*s 
and I's of the function to which it is ^plied without altering 
their number. 


4.1.2 A Special NRQF ; 

A special NRQF (s-NRQF) j^Karpovsky, I976j amongst the 
set of p-HRQFs of n =? 2t variables iss 


f®(x) 


"^i^i+t 


(4,2) 




.© x^X2t 


f (x) can also be represented in terms of the t x t Paley- 


Matrix P# 

' « ' 
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(x) = ig (1 -p ) 

^T^,e(x) 

where x = (x 2 ^ , x^ x^) , 

r (x)' “ x^/ . . . . /X^) is the bit- reversed form of r (x) / 

the right-half of xj 
r (x) = (x^# * * * ^ ' * 

e(x) = ' ^2t-l ' ••*•' ^t+l^ lef-half of x, 

P , is the (r (x) y e(x))-th entry of P. 

r (x) , e (x) “ 

Example 4.5 ; For n = 4 variables/ 
f^ (x) = X^X^ ®3C2X^ 

- % ( 1 - p ^ (x^X^) ^ 

s 

f (x) calculated using both these definitions is presented 
in Table 4»3« 

All Other NRQPs can be obtained from the s-NRQP using 
the transformation of Lemma 4.3* 

4,. 1. 3 Elatentary Properties of NRQFs 

PROPERTY i t f ® (x) ' =s Q for. Oi^ X:^ 2 ^, n = 2t variables, 

PROPERTY 2 i If f (x) is an NRQP of n = 2s variables. Then/ 
the number of ■ pnes' ^'in . f^Cx) , =■ 2^®"^ - -2®“^ and 

2s— 1 S** 1. 

number of zeros s= 2 + 2 

PROPERTY 3 1 Let N (n ss 2s) be the number of p-NRQFs of n variables. 
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Table 4.3 


SPECIAL NRQP OP 4 VARIABLES 
f® (x) = x^x^ ®X2X^ 


X 

(x 

4^ 

3^ 

2=^4 

x^.x^ 

X2.X4 

(x) 

r (!x) 

e (x) 

%(1-P. ) 

r (x) / e (x) 

0 

0 

0 

0 

0 

0 • 

0 

0 

00=0 

00=0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

10=2 

00=0 

0 

2 

0 

0 

1 

0 

0 

0 

0 

01=1 

00=0 

0 

3 

0 

0 

1 

1 

0 

0 

0 

11=3 

00=0 

0 

4 

0 

1 

0 

0 

0 

0 

0 

00=0 

01;=1 

0 

5 

0 

1 

0 

1 

1 

0 

1 

10=2 

01=1 

1 

6 

0 

1 

1 

0 

0 

0 

0 

01=1 

01=1 

0 

7 

0 

1 

1 

1 

1 

0 

1 

11=3 

01=1 

1 

8 

1 

0 

0 

0 

0 

0 

0 

00=0 

10 = 2 

0 

9 

1 

0 

0 

1 

0 

0 

0 

10 = 2 

10=2 

0 

10 

1 

0 

1 

0 

0 

1 

1 

01=1 

10 = 2 

1 

11 

1 

0 

1 

1 

0 

1 

1 

11=3 

10=2 

1 

12 

1 

1 

0 

0 

0 

0 

0 

00=0 

11=3 

0 

13 

1 

1 

0 

1 

1 

0 

1 

10=2 

11=3 

• 1 

14 

1 

1 

1 

0 

0 

1 

1 

01=1 

11=3 

1 

15 

1 

1 

1 

1 

1 

1 

0 

11=3 

11=3 

0 



■'4x 


0 

2 

'3 


0 X 2 

11 1 

1 1 -J* 

1 ^1 1 

‘■tmt 1 


3 

1 

1 

1 

1 
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Then N( 2s) = (2s-l) .N(2s-2) 

= (2s - l)i 

Each p~NRQP generates (2^®-l) np-NRQFs. 
Thus, total number of NRQFs of 2s variables 

= (2)^-^^ (2s-l)i 

(s~l)i 


<4»3) 


(4„4) 


4.1.4 Generation of NRQFs 

An NRQF Of n = 2s variables is generated by a combinatorial 
circuit of ANDs and EXORs getting its inputs from a mod- 2^ 
counter. 

Number of AND gates required = s. 

Number of EXOR gates needed = s-1. 

Example 4,6 : A typical circuit to generate 

f(x) a= x^X 2 ®X 3 X^®XgXg is shown in Fig. 4.1. 



FIG. 4,1 
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Conplexity of an N-length sequence vis-a-vis maximal— length 
sequence is guaged in terms of the number of stages of a minimal 
LFSR cxxcuit needed to generate it. This follows from the fact 
that only 2r terms of an r-stage LFSR- generated m-sequence are 
sufficient to know it completely. For an arbitrary sequence, 
longer is the minimal LFSR, greater is the number of terms 
necessary to decipher it completely and hence, more complex 
is it with respect to the same length m-sequence. Following 
this definition of complexity for NRQFs also, number of stages 
(r) of minimal LFSR to generate p-NRQFs of 2, 4, 6 and 8 varia- 
bles have been computed. Massey's shift-register synthesis 
algorithm [kassey, 1969J is used. The p-NRQFs and their 
respective r's are presented in Table 4.4. The notation 
adopted to represent p-NRQFs is that only their indices are 
given e.g. f(x) = x^X2®X23c^®X5X^®XgXg is written as 12 34 5? 68* 

T abl e 4 « 4 

■ n p-NRQF r 

' 2 -; ^ ■ 


^ p-NRQF X 

_ 




'■ 12 35 46; 

37 



13 24 56 

38 



13 26 45 

30 



14 25 36 

36 



15 34 26 

34 


' .34 ■ 

15 36 24 

37 


'34:, 

16 35' 24 

34 

' 32 .45 B 

34 ' 
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E 

)-NRQF 


r !1 n 

II . 

E 

)-NRQF 


r 

12 

34 

56 

78 • 

144 8 

12- 

34 

57 

68 

148 

12 

34 

58 

76 ' 

132 

12 

35 

46 

78 

15 2 

12 

35 

'47 

78 

140 

12 

35 

48 

67 

124 

12 

36 

54 

78 

156 

12 

36 

57 

48 

133 

12 

36 

58 

47 

141 

12 

37 

56 

48 

133 

12 

37 

54 

68 

136 

12 

37 

58 

46 

141 

12 

38 

56 

74 

129 

12 

38 

57 

64 

129 

12 

38 

54 

67 

129 

13 

24 

56 

78 

144 

13 

24 

57 

68 

150 

13 

24 

58 

76 

134 

13 

25 

46 

78 

152 

13 

25 

47 

68 

142 

13 

25 

48 

67 

126 

13 

26 

54 

78 

158 

13 

26 

57 

48 

131 

13 

26 

58 

47 

139 

13 

27 

56 

48 

131 

13 

27 

54 

68 

136 

13 

27 

58 

46 

139 

13 

28 

56 

74 

131 

13 

28 

57 

64 

131 

13 

28 

54 

67 

131 

14 

32 

56 

78 

144 

14 

32 

57 

68 

151 

14 

32 

58 

76 

135 

14 

35 

26 

78 

158 

14 

35 

27 

78 

140 

14 

35 

28 

67 

131 

14 

36 

52 

78 

156 

14 

36 

57 

28 

131 

14 

36 

58 

27 

135 

14 

37 

56 

28 

131 

14 

37 

52 

68 

142 

14 

37 

58 

26 

135 

14 

38 

56 

72 

132 

14 

38 

57 

62 

133 

14 

38 

52 

67 

133 

15 

34 

26 

’'^8 

158 

' 15. 

34 

27 

68 

143 

15 

34 

28 

76 

131 

15 

32 

46 

78 /■ 

152 

15 

32 

47 

68 

143 

15 

32.' 

48 

6 . 7 :''' 

127 

15 

36 

24 

78 

156 

; .15 

36: 

27 

48 

137 

15 

36 

28 

47 

131 

.;15..' 

;3.7' 

26 

48 

, ■"136 

15 

37 

24 

68 

143 

: : :15:^ 

37'' 

28 

46 

4:,/’ 131 

15 

38 

26 

74 

133 

15 

38 

27 


132 

15 

38 

24 

67 

133 

16 

34 

MM-' 


159 

16 

34 

57 

28 

131 


contd 
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p-NRQP r 51 n p-NRQF 

» - - ..* 


16 

34 

58 

72 

132 

8 

16 

35 

42 

78 

159 

16 

35 

47 

28 

131 


16 

35 

48 

27 

137 

16 

32 

54 

78 

159 


16 

32 

57 

48 

130 

16 

32 

58 

47 

138 


16 

37 

52 

43 

136 

16 

37 

54 

28 

131 


16 

37 

58 

42 

134 

16 

38 

52 

74 

133 


16 

38 

57 

24 

133 

16 

38 

54 

27 

132 


17 

34 

56 

28 

130 

17 

34 

52 

68 

142 


17 

34 

58 

26 

132 

17 

35 

46 

28 

130 


17 

35 

42 

68 

140 

17 

35 

48 

62 

137 


17 

36 

54 

28 

130 

17 

36 

52 

48 

137 


17 

36 

58 

42 

134 

17 

32 

56 

48 

130 


17 

32 

54 

68 

136 

17 

32 

58 

46 

138 


17 

38 

56 

24 

133 

17 

38 

52 

64 

133 


17 

38 

54 

62 

133 

18 

34 

56 

72 

130 


18 

34 

57 

62 

130 

18 

34 

52 

76 

130 


18 

35 

46 

72 

130 

18 

35 

47 

62 

130 


18 

35 

42 

67 

130 

18 

36 

54 

72 

130 


18 

36 

57 

42 

130 

18 

36 

52 

47 

130 


18 

37 

56 

42 

130 

18 

37 

54 

62 

130 


18 

37 

52 

46 

130 

18 

32 

56 

74 

130 


18 

32 

57 

64 

130 

18 

32 

54 

67 

130 








As exfjeicted from the non-linearity of NRQPs/ they are 
considerably more complex than the same length m-seqaences. 

A rough estimate of the complexity ratio from the given data is 
that the ratio is greater than 2^/2 .. _ 2 ’^""^ 

l.m Qm f • 

n n 

With these preliminaries on NRQFs/ we pass on to finding 
their WHT. 

4.2 WALSH- FIADAMARD TRANSFORM OP NRQPs 

This part is devoted to finding the WHT of NRQFs in 
general and will be useful in their correlation analysis. From 
now onwards, NRQFs will be taken in their converted 
form using the transformation: 

°f(x) = 

For the WHT analysis# only p— NRQFs have been considered 
because the WHTs of np-NRQFs follow directly from their 
respective p-NRQFs as follows; 

Let f{k) be a p-NRQF of n variables and g(k) be the 
np-NRQF which i a an sth dyadic shift of f (k); 
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<i,k> = *••• 


n 


^ i.k 


t=l 


t t 


is the ith LBF as pointed out in the last section. 
Substituting for'^i,k)> and f(k)^ their polynomial 


representations/ we get 


n 


^i k ® k k 

S (-1)^=1 ^ 


N~1 


k=0 


m, j 


m, 2 , nj- 


It is inferred that is the difference between the 
number of zeros and ones in the composite function (k) formed 
by the mod 2 sum of the ith BBF and the NRQF. 


Ci (k) 





9 S 


m J 


So, to know 


the number of ones and zeros in Cj^ (k) should 


be known. Each is formed out of an LBF from the linear 

vector space of LI^Fs of, n variables as i varies from 0 to N~l. 
Ihree repres^tative cases , of i are considered and the rest are 
shown to be combinations of these. 


■Case ,1'*^ : ,.i=0 ■ /f'' 

Case'''.2:'v‘Oh3..f Is one and the rest are ■ 


■Case 'Cl: 


.the ± 4 .*s 


t. 


5= 1, 2, . . , n) / say i^ and i^ are ones 
gp^iadiratic term of the NRQF f(k). 



Case 1: 


i =0 


P 


o 


N-1 

E 

k=0 


(- 1 ) 


f (k) 


= (number of zeros in f (k) ) ~ (number of ones in f(k)) 

_ ( 22 s -1 4 . 2 ®”^) - , (from PROP. 2 of 

NRQFs) 


SO/ 



■(4.5) 


Case 2 ; Only one i^ , for t = l, 2,.-.,n is one and the rest 
are zeros, 

n n-l p +1 p p -1 1 

Say, i=(0 0 ...0 1 0.. .0 ) , i.e. only i = 1. 

' ir 

or, i = 2 ^“’^. 

so, C. (k) = ^ i.k. © ^ k k . 

t=l m;tj 

m, j 




m, J 


The NRQP f{k) contains exactly one quadratic term, say 

•fm . 

cont^nirjg 


Thus, (k) = 



k k. 
w J 



68 

This function is nothing but a dyadic shift of f(k), say 
f(k®a) where 

^n ^n-l***\+l ^q-l“*^l 

a ={ 0 0 ...0 1 0 ...0) 

= 2^-1 

So, (k) = f (k® a) 

Thus#C^(k) is also an NRQP, 

So, for this case too 

= 2® ; i =2^"^, p=l, 2,...,n (4.6) 

Similar is the case when more than one i^=l / t=l, 2, ...,n, 
and each k^ for which i^= 1 belongs to a different quadratic term 
of the NRQP f (k) , i.e. at most one component k^ of each quadratic 
term forms the ith LBP. 

If t i, i. are ones and the rest zeros, where 

1 ^2 S 

t^, m , r tj^(^ 1, 2, • * • , then k^ , k_j_ , • • . should belong 

1 2 b 

to different quadratic terms. If such is the case, then each 
tern of the LBP can be combined with the ^propriate quadratic 
terms as shown before and the resulting (h) is another dyadic 
shift of f (k) * 

Example 4.7 ; Consider the NRQP 

f(k) = ^1^2® ^3^4® ^5^6 of variables. 

' ' ■ ":lo lo j o i> 

• - ■ ^ 5 ® 
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= k.® k.k-© k-,® k^k,® k^© k..k^ 

I 1 6 34 5 56 

= kj-k^ 

1 2 3 4 5 6 

= f (k ® (lOlOlO)) 

= f (k ©42) , 

which is the 42nd dyadic shift of f(k). 

Hence, "^21 

Thus, =2^ for i's of the above type. 

Number of such i's can be found out. There are s quadratic 
terms. So, number of such i's 
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(a) i's having one 

(0001) = 1 (0100) = 4 

( 0010 ) = 2 ( 1000 ) = 8 . ■ 

(b) .1 ' s having -two ' 1 ' ' s ; 

(0101) = 5 (0110) = 6 

(1001) =9 (1010) * 10 

So, for all these i's, = 2® = 4. 

These i's can also be enumerated using a graphical 
representation of an NRQF. Let n vertices of a planar graph, G^, 
represent the n variables of NRQF f(lc). v is sn edge of 

, the graph of f(k), iff Kpkg is a term of f(k). Hence, ciny 
NRQF can be graphed onto a disconnected graph with n vertices 
and s =n/2 edges, each representing a quadratic term. 

Example 4.9 ; Let f (k) = k^k^ ® ^2^4 

is shown in Fig, 4.2(a) 

£ 



n =4, No. of edges » 2, _ 

Corrplement of defined now: has all 

possible edges v^^v^ such that v^v^ is not an edge of G^, 
a# e € “s 1# nl-* 
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Gf U the complete graph of n vertices. and U 

for the f (k) of the previous exarrple are shown in Figs. 4.2(b) 
and 4.2(c), respectively. 



c 


G 


f 


PIG. 4.2(b) 



A path in a graph G is defined as a sequence of alternating 
vertices and edges with no vertices repeated. Length of a 
path with (L+1) vertices is L. 


Example 4. 10; In the '^G^ of Pig. 4.2(b), a 3-length path is: 

jTv^, ^ 2 ' ^2^3^ ^3^ "^^3^4' ^4j 

L 

It can also be represented as 


To enumerate all the i's of case 2, find all L-length 
paths (L «0, 1, 2, ,..,s-l)of ^G^ which are not permutations of 
each other and do not include vertices which are adjacent in 
G^. Each of ' theS'C ' paths ''represents ^ 'Sn ^ which F^=:2 in 

accordance with case 2. This will be illustrated with an 


4,11 : Let f(k) » k.k ® k 3 k^ ® k^k^, s = 3 


G^ and are shown' in Pigs. 4.2 (d) and 4, 2 (e), respectively. 
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PIG. 4.2(d) FIG. 4.2(e) 


For L =0, 1/ the L-length paths and corresponding i's are 
given in Table 4.5. 

Table 4.5 

0-length paths (igi^i^i2i2ij_) i I 2-lengths paths (igigi^i2i2^2^ 



0 0 

0 

0 

0 

1 

1 


kj. 


0 

1 

0 

1 

0 

1 


0 0 

0 

0 

1 

0 

2 


kj. 


0 

1 

0 

1 

1 

0 

. l>'3} 

0 0 

0 

1 

0 

0 

4 

{^1' 

^4.' 

’^s} 

0 

1 

1 

0 

0 

1 


0 0 

1 

0 

0 

0 

8 

{^2' 


’'sj ■ 

0 

1 

1 

0 

1 

0 


0 1 

0 

0 

0 

0 

16 

{ki. 

^3' 


1 

0 

0 

.1 

0 

1 

{>'61 

1 0 

0 

0 

0 

0 

32 

[^ 2 ' 

k3. 

’'el 

1 

0 

0 

1 

1 

0 

(These are all 

vertices 

of 


h^i' 



1 

0 

1 

0 

0 

1 

% >• 

i: 








k^/ 

>^a] 

1 

0 

1 

0 

1 

0 

1- 1 en g^H^p aiSi “ 







l-Imgth paths 







- 3 } 

0 0 

0 

1 

:0 

■1: 

v.,5 

1^3' 

^5} 

0 

1 

0 

1 

0 

0 : 

{^2' ^3} 

0 0 

0 

1 

1 

0 


{> 


kj} 

0 

1 

1 

0 

0 

0 

- ^4! 

0 0 


0 

0 

1 

9 

1= 

"1' 

’'sl ■ 

1 

0 

0 

0 

0 

1 

{v ^ 5 }^/: 

0 1 


0 

0 

1 

17 

r)c3. 


1 

0 

■Q:. 

1 

0 

Q 

(^2' 5^4}',/; 

iSiK 


% 

M 




^2" 


1 

0 

'O'; 

0 

1 

0 

4^2"' ^sl ' 

lUir 

11 



0 

' 'k- 

■rAB ■ 

i' 

^4' 

^s] 

1 

0 

1 

0 

0 

0 


These are edges of *^0 
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If path P, otherwise 0. 

N. = n + E ( ^ ). 2 ^ 

j = 2 J 

= 6 + ( 2 ) . 2 ^ + ( 2 ) . 2 ^ 

= 6 + 12 + 8 = 26. 


There are 6 zero-length paths (these are just the vertices of 
G^) , 12 one-length paths (edges of G^) and 8 two-length paths, 
Hence, the tabulated i's are the case 2 i's for f(k). 


Case 3: i is such that only two of the i's (t=l, 2, . . . , n) , say 
ip and i^, are ones and is a quadratic term of the NRQF 

f(k) 


^n ^n-l***^p4-l ^p ^p-l***^q+l ^q ^q-1*** ^1 
Let i=(00 ...0 10 ...0 1 0 ... 0) 


So (k) » k^ ® k^ ® k^k„ ® £ k,k 

m;^j 

m, jj^p,q 


^ -p - -q - "p"q ' — -j-m 


sa Kpk^ ® k 


i 1 ® ^ k .k 


£ ^ ./ kg, k ^ ) 


n' n-*r 


«:kg' . 3 ^,: :# 


j m 


m/ j/p, q 






a a , . a a . . 
n p+l p p-1 

;o * . . 0 , 1 0 


'• H _ • • « • '3. 


q 


. 1 


1 

. .0 ) 


This function is the conplement of the ath-dyadic shift of 

nuinbep'h*^ ones and zeros get interchanged vis-a-vis 
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so, = -2= (4.9) 

Example 4.12 ; Let f (k) = k^k^ ® k 3 k^ , n =4, 3=2'. 

For i = 12 (1100), 

= ^^^^2 ©k^k.^ ® 1 

so, F^2= -22 = - 4 . 


Cases 3 and 2 can be generalised further to two alter- 
natives: 

(a) i is such that the components of an odd niomber of quadratic 
terms and single components from any of the other quadratic terms 
form the LBF. This is similar to (4.9) and hence F^ = -2^. 


Example 4. 13: Let f (k) = ^1^2 ®^5^6 ®^7^8 ^ 

Consider i » (lljltllOl) = 253. 

So, C 253 (k) 5= kj© kg© k^® k^® kg© k,^® kg® f(k) 

sa k.® kjk,^ ® . kj ® kgk^ ® kgkg ® k,ykg 

■jss3' J ■ 




case 3 


a* , is -© ,;,IC © ^'5-^g © ^7^8 © ^ 

= f (k® (11111110))® 1 


-as':'" 


i is such" that the c 


cci^^nQits of 'an even number of 
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quadratic terms of f(lc) and single components from any of the 

other terms form the LBF, then P.= 2®. 

1 

g.g,. a S0Bl e .,4,,14; Let f (k) = k^k 3 © k^k^ © k^kg , n=6. s=3 
Consider i = (llllio) = 62 

“ ^2 ® ^3 ® ® ^5 ® kg © f(k) 

= kg © k^kg © kg © k^ © kg © kg © kgk^ © kgkg 

case 2 case 3 

= k.k„ © k^k, © Rck^ 

1 3 2 4 5 6 

= f (k ©59) 

So, Fggs 2^ = 8. 


Lemma 4,5 ; 

The unnormalised WHT of an N~point (N=2^®) NRQF, f (k) , 


is: 




, i s* 0 
, otherwise 


(4. 10) 


Bend©, in as''EKuch as each F^ ■ is -the distance (distance being 
the correlation) between the ith LBF and the NRQF f (k), every 
NRQF is equidistant from each elemaat of the set of LBPs, 


As an illustration for all the cases# consider an 

ggannle 4, IS; Let'.f'<k) "s* ® ^ 3^4 
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The u-WHT of f(k) alongwith the cases for each i are 
given in Table 4.6: 


Table 4. 6 


i 

(i 

4^ 

3^ 

2^1^ 

^1 

Case 

0 

0 

0 

0 

0 

4 

1 


0 

0 

0 

1 

4 

2 - 

2 ' 

0 

0 

1 

0 

4 

2 - 

3 

0 

0 

1 

1 

-4 

3(a) 

4 . 

0 

1 

0 

0 

4 

2 

5 ^ 

0 

1 

0 

1 

4 

2 

6 

0 

1 

1 

0 

4 

2 

7 

0 

1 

1 

1 

-4 

3(a) 

8 

1 

0 

0 

0 

4 

2 

9 v-- 

1 

0 

0 

1 

4 

2 

o 

c 

1 

0 

1 

0 

4 

2 

11 

1 

0 

1 

1 

-4 

3(a) 

12 

1 

1 

0 

0 

-4 

3 (a) 

13 

1 

1 

0 

1 

-4 

3 (a) 

14 

1 

1 

1 

0 

-4 

3 (a) 

15 

1 

1 

1 

1 

4 

3(b) 


Functions, which are conplements of NRQFs, also have 
similar u»WHT, 


If 'iCk) 

mm °f (k) 


k=0 ^ 
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N-1 
2 (-1) h 
k=0 


ik 


f(k) 



Hence, G = -F 

jLiSgS,§.._4.6; u-WHT of complement of an NRQP, f(k), is the 
negative of the u-WHT> F, of f(k). 

An important fact which stands out is that like the FZC 
sequences of previous chapter which have uniform-magnitude 
Fourier spectrum, NRQFs have similar Walsh— Had am ard spectrum. 
Accordingly, it is to be ejspected that NRQFs have ideal dyadic 
auto- correlation as the FZC sequences have ideal cyclic auto- 
correlation. This is demonstrated in the next section. 


4.3 DYADIC AUTOCORRELATION OF NRQFs 

Dyadic auto-ccrrelation of an N = 2^-point (n = 2s) NRQF 
f(x), f(x) being taken as *^f(x), is: 


N-1 

b-(s) <ss 2 




or fej = i I (2 .£) , 


)SOi 


' £'■ ''%;/4e :' ,thte':a-WjHT 'o£ ■ ’^f (x) 


/^s I 2^ J ^ 
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Thus, ^ H . (N N)^ 

~ N • ^ ^ Sq ' is the oth column of H . 

= (N 0 0)^, (H. H.^= N S . .) ■' 

-1 -j "ij'' 

This is the ideal dyadic auto- correlation. Hence NRQFs 
are ideal for LDSIV system identification. It is evident that 
any function having 'uniform.— magnitude WHI spectrum has ideal 
dyadic auto— correlation. This fact will be further dwelled 
upon when such functions called bent functions are presented 
in the next chapter. Presently, ,we go on to the cross— correla- 
tion properties of NRQFs. 

4.4 DYADIC CROSS-CORRELATION OF NRQFs: 

Dyadic cross-correlation between two n(=2s) variable, 

N (=c2^) -point , NgQFs f (x) and g(x) is defined as: 

® g(x®t) 

^ X«0 

i.e. the 't-shift dyadic cross-correlation is the difference in 
the number of aeros and ones of the coirposite function, 
f (x) ® g (x ® t) - the iitteer combination of\ f (x) and ' t - dyadic 
>hi fted 'fix)' * \ ^ ^ 

It is evident that# 

Let d(x) «= f (x®a) be tho resulting from ath dyadic- shift 

of its p-NRQF f(^x^,gCxJ foe another p-NROF* 



N-1 

Then/ b, (t) = (_2) ^ 

^ x=0 

Let x®a = k 

N-l 

So, b, (t) = 2 (- 1 ) f 

^ k=0 

= b^g(a®t) 

This can be stated as: 

Lemma 4,7: 

Dyadic correlation between ath dyadic shift of a function 
f (x) axid another function g(x) is ath dyadic shift of the 
dyadic correlation function b- (t) . 

This is nothing but dyadic shift invariant property of 
the correlation operator. 

Hence, we need to consider only the principal NRQPs and 
results for non-principal MRQPs follow in accordance with 
Lemma 4.7, , , ■ 

In general, two p— NRQPs, f (x) and g(x) , aare related 
thrO'Ugh pe.rmutatibn , of variolas as shown in Section 4.1.1. 
Permutations aro represented by their cycle structures.. As 
each yariable occurs exactly once in each of the NRQPs, their' ' 
cycle structure can be established just' by inspectibn. 

Example 4. ist Let ■' 
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be two p-NRQPs of 6 variables. 

(The notation shown in parenthesis is adopted to represent 
NRQFs) . 

The cycle structure ; 

(l) (264) (3) (5) 

consists of 3 single- length cycles or identities and 1 three- 
length cycle. This is called a one 3-length cycle permutation. 

Similarly, cycle structure of any two p-NRQPs can be 
established. We will base our enquiry into cross- correlation 
of NRQPs on the types Of cycle- structures relating them. Firstly, 
we will find cross- correlation between two NRQPs related through 
one 2-laigth cycle. Such permutations are also called Trans- 
positions. All the rest follow easily from these. 

4,4,1 Cross-correlation Betwea:i Transposition-Related p-NRQPs 
Let f (x) and g (x) be two N=2^-point (n = 2s) p-NRQPs 
related through a transposition. Such NRQFs differ in exactly 
two of their s quadratic ''tems. 

Exam &l^ ^ Let f (x) 

g(x) = x^X 2 ®X 2 X^®X 5 Xg 

Permutation cycle structure is (l) (23) (4) (5) (6) which is a 
.^Transposition, ; 

They differ only in ^their first two quadratic 'terms. 
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From the N-point functions^ f(x) and g(x)# construct 
the function r(x/t): 


r (x, t) = f (x) ® g(x ® t) 


So, bfg (t) = (number of zeros in r(x, t)) 
(number of ones in r(x, t)). 


In general, g (x®t) can take three different forms according 
to the value of t as has already been sho'^mi in Section 4.2. 

(i) g(x®t) = g(x) , t =0 (case 1 of Sec. 4.2) 

n 

(ii) g(x©t) =g(x) ® ^ c^x^ , (case 2 and 3(b) of Sec. 4.2) 

m=l 

c = depends on t and g (x) . 

n 

(iii) g(x®t) = g(x) ©"§1 c x. ® 1 (case 3(a) of Sec. 4.2). 

m=l 

Example 4.18 ; Consider g(x) = x^X 2 © ^ 3^4 
(i) Let t = (0001) 

So, g(x ®t) = ’^1^2 ®^3^4 

=s g(x) ® ^2 

■ ■ 

Ihus, cs=( 0010 ) tod = Xg (case 2 ) . 

I#! ' • 

' (ii) ' ■ 
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Thus c = (0011) , 


m 


m m 


Xi$X2 


(iii) For t = (nil) , 


g(x ® t) = ® X3X4 


g (x) ® X^®X2®X3®X4 (case 3 (b) ) 


Thus, c = (nil), g c^x = X ©x^®x^®x, 

m 


^2 ""3 ""4 


Hence, r(x, t) can take three different forms which are 
treated separately below. 

(A) r(x#0) = f (x) ® g(x) 

(B) r(x, t) = f(x) ® g(x) c^x^ 

(C) r(x, t) = f(x) ® g(x) ® 1 


Case A : 

Let the combination of the two quadratic terms in which 
f(x) and g(x) differ be; 

i.e. , the transposition la (l) (2) . . . (Jk) . . . (n) 

-rCsnO) ^ : ' v ' 
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Let r. (x) = X, ® x„ and 


2 ^^' ~ ® be LBFs of X/ 


(x) = 

X — { ^n, * • • 

assuming n^ i^ q'^ k ^ 1. 


So, r (x, 0) = r^ (x) r^ (x) 


N-1 


Number of ones in r(x,0) = s r(x,0) 

X=:0 


b = 

= N H (R, . R.) 

, where, ' 

“"IT^ IT ^ 
12 

ss -^1 ' — 2 


) (0) = 2^® 

N-1 

Thus, fc 

2 (B 



m=0 

So, 

and R^ have to 
~2 ■ 

be fotind 

WifT, J, 

of the/cth N*point LBR 


N-1 

S r (x) r_ (x) 
x=0 ^ ^ 


= b (O), the zeroth dyadic 

correlation between 
r^ (x) and r^ (x) 


-2 N i -2 * 


X -2 m 


m 


m m 


y ea- 
rn 


% ; yVm 






‘vl?: 




fitm, ij ^ ' f • 

' , - ' . . ' ^ ':'"V/"' ' 'i/V ■: : ^'■: ' 
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and 


( 0 


n 


^c. Ic 


Similarly, ^2^^^ 


... 1 0 . 

n 

= 

m=l 

= Xj * 


q 

1 0 .. 


1 

. 0 ) 




(4. 12) 



2 

2 



2 

2 

®n 




®k 

c = ( 

0 . . 

. 1 

0 

m m m 

1 0. . 

0 

1 


1 

c 


N-1 

So, R^ = ( ^5 

0 ... 

0 

-h 

0 

. . . 0 ) 

0 

1 


2 

c 


N-1 

and R 2 ( h 

0 . . . 

0 

-h 

0 

. .. 0 ) 


0 

1 



N-1 j 

Thus, H^.R^ = 

( 1/4 

0 . 

. . . 

• * * 

.. 0 ), 


1 / 2 


(4. 13) 


N-1 N-1 

So, 2 r(x,0) = N S (R. .R_)^ 

x=0 mcrO ^ ^ 

= 22 ^ 1/4 = 


O a O 

Or , number of ones of r(x#0) = 2“^ “ . 

Number of zeros of rCx, 0> > = 3.2^®“^ 

Di f f erenoe ■ » ^ ^ • '= . ' 

Hence, b^ (O) * 2^®*"^ (4,14) 

i.e*,, ■ Idle: trmsposition- 
related 2s-varl5ibie .^NRSliE'S' is 2 * 
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§ • Before we go on to find b^g(t) for this case/ we note 

the relationship between ^c, ^c, c and t. 

(i) If t is such that c = ^c, then 

g(x e t) = g(x) ® S^c x^ / ^c is as in (4.12) 

m 


= g(x) & X. 9 X 
= g(x^...Xi-..Xq..,Xj 




n 


SO/ t = ( 0 ,, 


2 

= c 


1 0 


1 0 . 


(from 4. 13) 


Xj,. . .x'^) / as 

XjXq are terms 
^ of g(x) - 

0 ) 


O 1 

Hence/ t = c <=> c=: c 


(4.15) 


(ii) If t is such that c = ^c, then 


g (x ® t) = g(x) 9 


m 


m m 


g(x) $ Xj ® 
g (x^« * • ^i * * * 


"'n : , . i , ' 1 

t '( 0 ... ... . % ' ;()«’. O'. . « .'* 0 ' ) 

ags ^ 'C, ' ' : 

1 ' 2 ' 

Thus, 


Xj) 


(4.16) 


.If:, "d 












.0 ) / 
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then g (x @ t) = g(x) ® X. ® X © x. ® 


X, 


g 


So, t = 


’ ' • * ^q* * * ' 

t t . t 

n "'1 

0 , 


■ . -x^) 


\ -^1 

. 1 0 1 0 . . . . 0 ) 


Thus t = 


3^ / V 3 

c <:=> c = c 


(4.17) 


For all the rest of t's, i.e. , for a t such that a^^c, or 

3 .A 2 3 

c, c, c or c. 

This will help in finding for this case and the 

next. 


For this case, 

r (x, t) = f(x) © g(x) < 


m 




(x. © Xq) (xj ® x^) 


® 


m 


m m 


= r^ (x) r^ (x) ®r^ (x) , r^ (x) = 


^ — m 

in 


S r^ 2 ®r2(x) , r^ 2 = r^(x)r2(x) 


N-1 


£ r(x#t) a £ (r^ 2 ^^^ ® r^(x)) 


x=0 


X 


* s (Arithmetic represen- 

^ tation of ffi) . 


» 'S 


©■ £ r^Cx) - 2 £ 


(4,18) 


From "ease A# 


£r,2(3c) a 2^®*“^^ 


X 


12 


(4,19) 



(x) is an LBF 


so E r_ (x) = 2 

X 3 


2s-l 


(4.20) 


£ r2(x)r3(x) = b „ (O) 

X ^ ^12^3 


(4. 21) 


To find this, we calculate R^2 and R^. 


(a) = 


(x) 




So, r . 


■^c 


(from 3 . 15 ) 


Similarly, r2 (x) = ^ '^c^ 35^ 

m=l 


So, £2 


(H - H. ) 
—O —2 

c 


-12 " -1^-2 


c c 

5 (Ho -«! - ^2 + H ) , h = h * 2^, H^.H =H 

c c c 


So, R 


± i %2 


' fi ' 4 '■ ^ ' ,-f ■' Ho ' ,) 

, 2 $ 4 ' 4 ': ■ ■■ : , ' - 3 ^,' 


C2s*2) [(j23; o__...o )'^ - 

0 N -1 
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(b) £3 


n 


r jCx) = ^ 


o 


N-l 


So, R 


3 


(1 0 -1 0 0 ) 


Novr, b ^ (0) = 2^® ,S (R_.R-). 

^12^^ 3 t=0 


N-l 


(4.23) 

(4.24) 


- 12-3 


h 


1 

0 

m 

♦ 

"1 

0 

* 

-1 

0 
« 
m 

1 
0 


0 




2 X 


4— 


1 

0 


0 

-1 

0 


4-0 


0 


(From (4.22) and (4-23)) 


^c, and '^c are fixed once f{x) and g(x) are known. As the 

dyadic shift, t, varies from 0 through N-l, it generates the 
set of all possible LBFs X %V ^ through 0 to N-l 

From, the nature of the product, £^2’-3 ' four distinct 

possibilities : for t 'is, 

■ 1 ' 

,, , {a^ cw, c ,, 

(b) cw,^c ' 

, Cc) c*^c 

■ , ■ ■ . (d) c/^c# oiT'.^e-’ ' 
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X§;1 ■ ■l . lsL * c = c — y t = ^Ct (from 4, 15) 

2 1 

c = c => t = c. (from 4.16) 

For both of these, 

0 1 c=^cor^c N-l 

3 ~ 8 ^ ^ 0 ... 1 0..,...0 ) 

So, from (4.24), 

b„ „ (0) = (1+1) = 2^®"^ 

^12^-3 

So, from (4.2l), 


( 4 . 25 ) 


Substituting (4.19), (4. 20) and (4.25) into (4.18), we get 


N-,1 


£ r (x,t) » 2 ^®”^ + 2 ^®“^ - 2 . 2 ^®“^ 


x*0 


ss 2, 


2 s** 2 


So, number of ones of r(x, t) ss 2 


2s- 2 

2s «2s-2 ,^2s-2 


number of aseros of r(x, t) «» 2 -2'"'^ ^ = 3.2 

Thus, the difference*: 3* 2^®" 

Hence, b^^Ct) » 2^®*"^, tat^c, 

t c ** c , ■ 


(4.26) 


So, t as C 




6: ■ i: ■&■ ri&kii'i , 


.',|L:,: 


T 
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12 3 


or 




Subscituting (4,. 19) , (4.20) and (4.27) into (4.18), 


we 


N~1 

Z t (x. 1-) 
x=0 


j2s.2 22s-l_g ^ 3 _j2s-2 


30. 

b.p, (g) 

EO , 


Ml 

"• C .'i s 


This 

So, 

” 12*^3 

or. 

b„ „ 
i7r3 


2^'*' - 7... 3. 2^®“^ 

-2S-1 


'2''^“"- , t =^c =^c ® 


1 2 


T 


(0) » 2 r^2^5c)r2(x) 

!}C 

2 3 

» 2 


N**** X '' ' ' 

SO', . r{x,t) * '2^^":'^: 




2s''**3 


x«0 


^ ' 2 ' '' 


b. 


at 0 








(4. 27) 

get, 

(4. 28) 


(4. .29) 
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some of the t's not equal to or will give rise 

to case C . 

n 

Case C ; r(x, t) = f(x) ® g(x) ® "W c x ® 1 

mm 

m=l 

= (x. ® X„) (x . ® X, ) © C X © 1 

^ <3 J K ‘-'T' m m 

m 

= ^ © 1 / (x) as 

before 


It is to be noted that this form cannot occur for 
12 3 

t = C/ c or c i.e. for this case c cannot be either one 

12 3 

of C/ c or c. 


As q(x) has quadratic terms x^^Xj, and ^jX^, for this case? 

if c » C/ then x^x^ has to be a term of g(x) which it is not? 
and if c « then XjX^ has to be a term of g(x) which it is not 
if c*^c. 


g(x) 


gCxj^- 


■ X, 


HI. # • . 


Xj 



« gCx) ® X£ © Xg © Xj © Xj^ 

which is not of; the of '' cas'e-G... So^ for this case, the 

only alternativa is of case B(d) i.e., t/ c, c or c. , 


for this as in case 

bji^Ct) *» 0 


(4, 30) 



from (4»l4)' 
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Lemma 4.8; 


Let f(x) and g(x) be two N-polnt (N=2^®=2^) p-NRQPs 
related through a transposition (l) (2) . . . (jk) . . . (n) . Let 
f(x) contain and g(x), ©x^x^. 


- n i q I 

Let c =( 0....0), 

2 n j k 1 

c™( 0 ****!l 0 »*** 0 ) and 

3 1 2 

c = 'C ® then. 


bjg(t) 


0 

L 


1 2 

t=0, c 


Otherwise 


Lemma 4.7 and 4*8 jointly lead to another interesting result. 


Corollary 4* It 

Cross- correlation between g(x) and ath dyadic shift 


of f (x), d(x) , Is 


Vi 

d©' 


Ct) 


2 , 2 ®- 1 ; ^ ^ ^ ® a, ® a 

'4'- 

j'r'.O ' ■ ,r ;,o€he,iwi se ; 


Hence, it^ fDl:|d^Sx:thilt^^;/■;;^X::;.^ 

; Cl) :nwnber 

(11) three of these are ; tod' -the fourth is -2^®“^. 

We illustrate this vrl'^ at exapt^le* 
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Exanipl^ 4 . 9 : Let f(x) = 

ana g(x) = XjXjSxjXjSx^x^ax^Xg 

bo two p-NRQFs with n = 2s =8, H = 2® = 256. The permutation 
cycle stxucture is 

(1) (2) (3) (4 5) (6) (7) (8), (j = 4 , k = 5 ) 

f(x) contains XjX^SxgXg (i=3, q=6) 
g(x) contains >£: 3 X 5 ®x^Xg . 

So, s (00100100) = 36 

a (00011000) » 24 
* (00111100) « 60 


Hence# 



# t *0# 24# 

f t * 60 

# otherwise 


36 


% 


Moving ■ found the , cross- oarr elation for transposition-related 
p-NRQFS'#, we will generalise these for arbitrary p-nrqFs. 


cro.s?-Coi^ye latf.<^; , ^^bjtr axry' P-NRQf^ 

■ , ' It: it: 

relating two fu«CtlonS./,.^'#;’‘ai5Bt';;g(^ _be obtained through '■ 
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a sorioG of transposition-roi 

related intermediate functions. The 

cross-correlations between each c-f 

cuwe^n each of these pairs can be found out 

using thu results of the previous section, some relation connect- 
ing those pairwise correlations is In order. This Important 

relation is derived below: 

beranTi_i2‘ Let fix), g(x) and r(x) be three NRQFs of n = 2s 
variables. Then, 


3X1 


In words#- the* cross-^correlstion between f and g# b is 1/N 

"Ly 

times the cross— correlation between the correlation functions 


~fr ~rg' 


N-l 


mOOFt R.H.S. E b.^(k) b_ (k®t). 

k»0 


In vector .notation, 


R.H.S. « ^ ^ H (Sfx-'*Srg^' 


where B * • H b- and 


S WfW ^ 'W. 'b'" m- ■ 

jETCI " w ” 


Now, m (£*R),andw o 

^ CS ^ ^ ^wlwveo,>' 




' * 2 i ^ s ' 
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So, R.H.S. * -3^ H (P.R.R.G) 

Now, R = (2® + 2s ±2®)'^ 

Thus, R . R = 2^® (1 1 1)^ 

SO, RtiS = -3^ . N. H (P . G) 

= 1 g (H • G) 

* ^fg Q.E.D. 

If two p-NRQFs f(x) and g(x) are related through a series of 
transposition - related p-NRQFs as: 


f Cx) , , . . . . ^ ' g (x) ; 


. ...fa sure the transposition - related 

####•• 

, b* ^ „ are known by the results of 


f, , £2 

pairs, then 

X' X ^ 

the previous section, furthet, by using Lemna 4.8 repeatedly 


■^1^2 




’na-l 


f m 


with the pairs » £• ^f^f '* *' 

: 12 12 2 3 

fef fefrr ^ found out. 

In general, if two 11 variable functions f(x) and gCx) are 
related through an li-lengt^i permutation cycle ^ 


{i j^i 2*' 


they oan in turn tm 
positionst 


rii^^'^;;;4a':;3s‘Olations ^ aiiKsngst (iXl) . :trans- 
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(ii ij,^) 

{il ij) 

(i 2 '^3^ / i i 2* • • 1/ 2/ . . . /n.^ 

^^L~2 ^L-1^ 


ExaiTiple 4. 20 : 

Let f (X) = # X ® XjXg ® x^Xg and 

g(x) = ^ 1^7 ® ^2^4 ® ^3^6 ® ^5^8 


The cycle stxuctxire is (l) (2357) (4) (6) (8). This can be obtained 
through the following series of transpositions: 

(1) (27) (3) (4) (5) (6) (8) . (f f 

(1) (23) (4) (5) (6) (7) (8) , (f^ f2) 

(1) (2) (3 5) (4) (6) (7) (8) , (f2 g) 

The resulting series of functions is* 


f ( x ) ■ ^ 1^2 ^ ^ 3^4 ® ^ 5^6 ® ^^8 

fj(x) » « X3X^ <6 XgXg ® X2Xg 

£2(x) • XjX^ #; 

g(x) « x^Xj 9 X 2 X 4 ® ^CgXg ® XgXg. 


Firstly# ^ g' 


£^Jt 




r, (x) • x^ t Xg r- ^ . 

So, ^c*« (10000003,) ml29 (O^OOOOIO) 

*1’ '(110000 11) • ' 1^5 
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128 , t =0, 66, 129 


(t) =<-128 , t=195 


ff. 


I 


, otherwise 


— f -f 
12 


(x) = X 2 ® x^ 


so, c= (00000 110) = 6 

3 


r 2 (x) = x^ ® Xg 


, “c = ( 10001000 ) = 136 
c = ( 10001110 ) = 142 


So, , (t) 

12 


128 , t= 0, 6, 136 

•128 , t = l42 

0 , otheJTwise 


/ t- j • 


(x) = Xg © Xg , 


r2(x) =Xg © Xg 


So, ^c= ( 00010100 ) = 20 , = (10100000) = 160 


SO, 


c = 


b. „(t) = 


( 10110100 ) = 180 

128 , t =0, 20, 160 
-128 , t = 180 
0 , otherwise 


An important point to note here is that for each transposition- 
related pair, (f (x) , f j (x) ) , ,e. , the set of 

dyadic shifts for which b_ ^ is non- 2 ero, form an Abelian group 
with the group operation being © * This is evident from Lemma 4.7. 

I 12 

Also, every distinct tran^Osition results in distinct c, c 



ya 


and c vis~a— vis another. As seen above/ no two shifts except 
0 for which — f f — f g non-zerO/ are same. 


Now/ 


'ff. 


- 1-^2 

(t) = - £ b 


KT ^-F-F ^ (x©t) 

x=0 1 1^2 


This irtplies that b^^ is non-zero for all those t's for which, 

X ffi t = X/ 

12 

where, x and x are the dyadic shifts for which b_^ and 
”^1^2 non- zero, respectively. 

So, t = ^x ® ^x , for all ^x and ^x. 


and each b^^ (t) = + — ~ (2^®”^. 

= + 2^-2^ 


1 2 

' + • if b^£ ( x) and b^ ^ ( x) both are of the saitve sign and 

1 ^12 

if they are of opposite signs. This follows from the two proper- 
ties mentioned above. 


Using b^/r and b^ ^ , non-zero b^^ (t)sare tabulated 

“l ^1^2 ^2 

in Table 4,7, 

For ail other t's, b^^ (t) ='■•0. 

So. b r-j- has 4 tim es (l6) the number of non— zero values of its 

f — ff 


component arrays magnitude is halfed 

to 2^^*"^ obviously#' These non— zero— b^£ dyadic shifts, t, also 

form a group because they constitute the "set of all possible ® 

suras of the elements of similar groups of each of the arrays 


b/:^ and b^ ^ . 

-*1*2 
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Table 4.7 


Figures in parenthesis are (t) 

12 ^ 
t = -^x X 


2x 

lx 

0 

6 

136 

142 

0 

0 

6 

136 

142 


(64) 

(64) 

(64) 

(-64) 

66 

l6 

68 

192 

204 


(64) 

(64) 

(64) 

(-64) 

129 

129 

135 

9 

15 


(64) 

(64) 

(64) 

(-64) 

195 

195 

197 

75 

141 


(-64) 

(-64) 

(-64) 

(64) 

of +ve 

values 

li 

o 

No, of -ve 

values 


= 3.3+1 ; , = 3.1 + 3. 1 = 6 


Simil arly. 


b 


fg 


(t) 



N-l 

‘’f,g 

x=0 2 2^ 


Again, using similar arguments as for b^^ , 

2 

dyadic shifts, t, for which b^^ is non-ze^'o. 


we compute the 
in Table 4.8, 


Again, the number of non-zero values has increased four- 
fold to 64 and their' magnitude has decreased by half to 32. 


Number 

of +ve 

values 

= 10 

X 3 

+ 6 

X 1 

= 36 

Number 

of -ve 

values 

= 10 

X 1 

+ 3 

X 6 

= 28 


The t‘ s again form a group 
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Table 4.8 


Figures in parenthesis are the correlations 



2 

X 

1 

X 


0 


20 

160 

180 

0 

0 

(32) 

20 

(32) 

160 

(32) 

180 

(-32) 

6 

6 

(32) 

18 

(32) 

166 

(32) 

178 

(-32) 

136 

136 

(32) 

156 

(32) 

40 

(32) 

60 

(-32) 

142 

142 

(-32) 

154 

(-32) 

46 

(-32) 

58 

(32) 

66 

66 

(32) 

86 

(32) 

226 

(32) 

246 

(-32) 

68 

68 

(32) 

80 

(32) 

228 

(32) 

240 

(-32) 

192 

192 

(32) 

222 

(32) 

106 

(32) 

126 

(-32) 

204 

204 

(-32) 

216 

(-32) 

108 

(-32) 

120 

(32) 

129 

129 

(32) 

149 

(32) 

33 

(32) 

53 

(-32) 

135 

135 

(32) 

147 

(32) 

39 

(32) 

51 

(-32) 

9 

9 

(32) 

29 

(32) 

169 

(32) 

189 

(-32) 

15 

15 

(-32) 

27 

(-32) 

175 

(-32) 

187 

(32) 

195 

195 

(-32) 

215 

(-32) 

99 

(-32) 

119 

(32) 

197 

197 

(-32) 

209 

(-32) 

101 

(-32) 

113 

(32) 

75 

75 

(-32) 

95 

(-32) 

235 

(-32) 

255 

(32) 

141 

141 

(32) 

99 

(32) 

237 

(32) 

249 

(-32) 


These results are generalised below: 


Lemma 4.10 ; If two n = 2s variable NRQFs, f(x) and g (x) ^ are 
related through t transpositions, then; 

t 2 1 

(a) Number of non-zero values of = (4) =2 , 

^ ' 2S“* 1. ' ■ ' 

(b) Magnitude of tbe values = ih) *2 

= 
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■o 4- *1 4” *1 

(c) Number of positive values = 2*^ + 2^’*'^ 

(d) Number of negative values = 2'^ - 2^~ , and 

finally, (e) all the dyadic- shifts, resulting in non- zero 

bfg(a) form a group with ® as the group operation. 

An s-length permutation cycle is the longest possible 
between two NRQFs, f(x) and g(x), of n = 2s variables. So, least 
cross- correlation magnitude achievable according to (b) is 

^2s-(s-l), t=s-l here. 

= 2^”*"" = 2(n)’^ , where N =2^® is the length of NRQFs. 

So, lim 2(N)'^^ _ lim _2 

N 00 N N GO “ • 

Hence, asyrrptotically-ideal ' cross-correlation is achievable 
with NRQFs, 

Representative examples for b^^ for n = 8 variables and 

various permutation cycle structures - viz. one 2-length cycle 

(one transposition), one 3-1 eng th cycle (two transpositions), 

two 2-length cycles (two transpositions) and one 4-length 
cycle (l;ihXbe transpositions) - arc given in Tables 4.9, 4.10, 

4. 11 and 4. 12. 

This brings us to the end of dyadic correlation analysis 
of NRQFs. Cross-correlation between NRQFs related through 
arbitrary permutations has been found out. Enumeration of all 
possible L-length permutation cycles for ^ NRQF will complete 
the picture. 



Table 4.9 


Functions f and _g related thru one 2-length permutation cycle 
f (x) ^ = x., X2®x^x^®XgX^®XgXg 

Cycle Structure : (l) (2) (3) (4) (5) (67) (8) 



-128- 240 

0 otherwise 


Table 4. 10 

Functions f and 3 related thru one 3- length permutation cycle 
f (x) = x^Xj^XgX^SxgXg^x^Xg g (x) =X2X^®x^x^@XgXg®X2Xg 


Cycle Structure : (l) (237) (4) (5) (6) (8) 


bj^(lc) 

1 k's 

1 



64 

0 6 9 

66 68 

77 129 135 136 202 

-64 

15 75 14 2 

195 197 

204 

0 

otherwise 





Table 4.11 

Functions f and 2 

rel ated 

thru two 2-length permutation cycle, 

f (x) = 


^6®^^8 

g (x) = x^Xg^XgX^^XgXgSx^Xg 


Cycle structure : 

(1) (25) (3) (47) (6) (8) 

i 

k's 




64 0 18 33 72 90 105 132 150 165 255 

-64 51 123 183 204 222 237 

0 otherwise 
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m 

where each L."^ 1 and s L. <' 3 . 

X . X 

1=1 

Example 4,21 ; For n =8 variables, consider the p-NRQF 

f(x) = X. X„®XnX^®XcX^©Xr 7 Xo 
12 3 4 bo 78 

Number of p-NRQFs of 8 variables = N (8) = 105 

" ’* single 2-length permutation cycles of f = 2 . (l)i=12 

" " single 3-length peririutation cycles of f = ~C^.2^, (2)1=32 

” " single 4-length permutation cycles of f = 2^ (3) i = 48 

" " tv70 2-length permutation cycles of f = "^^ 2 * (l)i) 

= 6 X 2 = 12 

This exhausts all possible p-NRQFs of 8 variables. 

Some computed results on the periodic and aperiodic corre- 
lation parameters of p-NRQFs are presented now. 

4.6 SOME RESULTS ON THE PERIODIC AND APERIODIC CORRELATION 
PARAMETERS OP p-NRQFs 

Let C (t) and p (t) be the aperiodic and periodic cross- 
correlations between two p-NRQFs u and v of length N =2^. 

p (t) = c (t) ■+■ c (t-N), 0 <t:< N-i 

•'=^uv UV UV •Ni: 

This is also called even crOss-correlation whereas odd cross- 
correlation, is; 
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p (t) = C (t) -C (t-N) , 0Xt<N-l . 

■^uv uv uv ' ' ^ sr* 

p^(t) and ^^(t) are the even and odd autocorrelation functions. 
Certain correlation parameters of importance in SS-CDMA systems 
[pursfey and Roefs/ 19793 are: 

(1) Peak even cross-correlation^ M(u,v), 

M(u,v) ^ ' 

( 2 ) Peak odd cross-correlation, M(u,v), 

M(u,v) ^ max"^ ‘ ° t^N-lJ'; 

( 3 ) Peak out-of-phase even auto- correlation, M(u), 

M(u) » max|^ 

( 4 ) Peak out-of-phase odd auto- cor relation, M (u) , 

'Kin) max|^ (p^(t)| : l^t^N-lj^; 

( 5 ) L (u) =11^ t;|p,^(t){ =M(u) , 1:^ t^, ' where {1 a|} 

is cardinality of A. L (u) is the number of occurrences 

of M(u) inp, (t); 

u 

(6) Si delobe energy, 

N-l „ 

S (u) = S C (t) 
t=l 

The relevant parameters for p-NRQPs have been found in the 
following way: 

(i) Auto-optimal (AO) phase of the p-NRQF is found out. 



The phase j of a sequence u Is an auto-optimal phase if 

^ (<p-J u) = p^Q (u) axid if L(T-'u) is the minimum of L(T^'u) over 

all i for which M(T^u) = ^^_(u) . 

— -^AO — 

In other words, M(u) is found out for each of the N-phases 

of the NRQF u and number of times it occurs in each phase is also 

computed. Then, minimtim of M(u) ever all phases is computed. 

Say, M (T^u) is this minimum. But there can be many j's with 
/\ 

this M (u) . Sd, those j*s are chosen for which number of 

/V 

occurrences of M(u), L (u) , are a minima. Even then there can 
be more than one such phases. These are the auto-optimal phases. 

(ii) Amongst the auto-optimal phases, those are chosen for which 
u has the minimum sidelobe energy, S(u). Now', only one or a fevr 
phases are left. This/these is/ are the AO/LSE p-NRQF. 

(iii) (i) and (ii) are repeated for each p-NRQF of n variables. 

(iv) Peak odd, and even cross-correlations, M(u,v) and M(u,v), 
are computed for each pair of AO/LSE p-NRQFs of n variables. 

(v) Peak even out-of-phase auto- correlation of each p-NRQF is 
computed. 

It is to be noted that, 

M(u,v) = M(v,u) and M(u,v)= M(v, u) 

Correlation parameters for all p-NRQFs of 4 and 6 variables 
and for some of 8 variables have been computed following the 

above procedure# Results are reported in' Tables 4. .13 ,and 4# 14# 
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Table 4. 13 


AO/LSE Phases of p-NRQPs of 4,6- and 8 Variables 


n 

N 

P 

-NRQF (u) 

AO/LSE Phase (s) 

(u) 

L (u) 

S(u) 

4 

16 


12 

34 

6 

2 

■8 

24 




13 

24 

3 


8 

24 




14 

23 

4,7, 12, 15 

6 

2 

40 

6 

64 

12 

34 

56 

26 

10 

4 

752 



12 

35 

46 

12 

14 

4 

9 28 



12 

36 

45 

43 

20 

4 

1480 



13 

24 

56 

20 

16 

2 

1088 



13 

25 

46 

21,53 

14 

4 

1272 



13 

26 

45 

21 

12 

10 

1096 



14 

23 

56 

21 

12 

4 

992 



14 

25 

36 

41 

12 

4 

1112 



14 

26 

35 

11 

12 

2 

1152 



15 

34 

26 

12 

14 

2 

992 



15 

23 

46 

54 

14 

2 

1088 



15 

36 

24 

11 

12 

2 

1000 



16 

34 

25 

9,41 

12 

2 

840 



16 

32 

45 

20,52 

. 10 

6 

816 



16 

35 

24 

9,41 

12 

2 

1064 

8 

256 

12 

34 

56 78 

102 

5Q 

2 

20720 



12 

34 

57 68 

47, 

56 

4 

21824 



17 

24 

56 38 

167 

48 

6 

28688 



17 

24 

36 58 

215 

52 

2 

29456 



15 

37 

26 48 

80 

58 

2 

31744 


Table 4. 14 


Peak Correlation Parameters for AO/LSE p-NRQFs 

(a) n =4, N = 16 


p-NRQP 

12 34 

^ 13 24 

14 23 

t 

12 

34 

4 

8 

s 

13 

24 ' 

10 

4 

8 

14 

23 “ 

1 

10 

8 

/I 



(b) 

n =6, N = 64 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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2 

2 

2 

3 

3 

3 

4 

4 

4 

5 

5 

5 

6 

6 

6 




5 

} 

3 

3 

3 

2 

2 

2 

2 

2 

2 

3 

2 

3 

3 

3 

3 

P- 


! 

3 

4 

5 

6 


5 

6 

2 

5 

6 

4 

3 

6 

4 

2 

5 




5 

3 

3 

5 

4 

4 

5 

4 

4 

5 

3 

3 

2 

4 

2 

2 

4 

2 




1 

J 

6 

6 

5 

6 

■ 6 

5 

6 

6 

5 

6 

6 

4 

5 

5 

4 

12 

34 

56 

"T 

1 

» 

1 

16 

32 

32 

32 

16 

16 

32 

16 

16 

32 

20 

20 

3 2 

16 

16 

12 

35 

46 

I 

1 

1 

36 

16 

32 

16 

32 

16 

16 

24 

32 

16 

32 

24 

16 

20 

32 

12 

36 

45 

1 

1 

1 

26 

30 

16 

16 

16 

16 

32 

16 

16 

16 

32 

20 

20 

32 

20 

13 

24 

56 

1 

I 

1 

24 

20 

28 

20 

32 

32 

32 

20 

16 

20 

16 

32 

16 

16 

32 

13 

25 

46 

1 

I 

1 

18 

26' 

14 

30 

16 

32 

16 

32 

16 

24 

32 

16 

32 

20 

20 

13 

26 

45 

8 

t 

i 

18 

16 

18 

30 

32 

16 

16 

16 

32 

32 

16 

16 

16 

32 

16 

14 

23 

56 

f 

1 

I 

30 

18 

32 

30 

22 

22 

16 

32 

32 

20 

32 

16 

16 

32 

16 

14 

25 

36 

t 

1 

i 

18 

20 

16 

18 

20 

20 

20 

20 

32 

20 

20 

32 

32 

16 

24 

14 

26 

35 

1 

\ 

20 

30 

14 

18 

20 

20 

24 

20 

16 

32 

16 

16 

16 

20 

32 

15 

34 

26 

i 

1 

1 

20 

16 

14 

18 

24 

30 

30 

18 

34 

16 

32 

32 

32 

20 

16 

15 

23 

46 

t 

f 

1 

24 

26 

24 

20 

18 

22 

26 

20 

24 

24 

16 

32 

16 

32 

16 

15 

36 

24 

1 

1 

1 

20 

24 

20 

30 

22 

16 

18 

22 

16 

34 

18 

16 

16 

16 

32 

16 

34 

25 

1 

1 

t 

24 

20 

18 

20 

18 

22 

20 

22 

20 

26 

13 

18 

16 

32 

32 

16 

32 

45 

f 

1 

1 

16 

20 

22 

22 

20 

30 

30 

20 

20 

20 

24 

18 

26 

16 

32 

16 

35 

24 

I 

1 

1 

I 

t 

20 

30 

20 

22 

22 

22 

18 

22 

36 

20 

20 

28 

32 

22 

20 


. . . contd. 
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(c) n =8/ N = 256 



J 

1 

1 

1 

1 

1 

1 


1 

S 

! 

2 

2 

■ 2 

2 

2 


1 

t 

3 

3 

2 

2 

3 


f 

S 

'•'r 

4 

4 

‘•■1’ 

7 

p-NRQFs 

t 

5 

t 

5 

5 

5 

3 

2 


1 

f 

5 

6 

7 

6 

6 

6 


5 

C 

7 

6 

3 

c; 

kJ 

tr 


i 

I 

1 

8 

8 

8 

8 

8 


; 






12 34 56 78 

1 

1 

64 

128 

66 

40 

64 

12 34 57 68 

i 

i 

146 

64 

40 

48 

40 

17 24 56 38 

J 

i 

56 

4 2 

64 

128 

48 

17 24 36 58 

5 

52 

48 

80 

64 

4t‘'r 

15 37 26 48 

( 

5 

i 

1 

56 

46 

48 

4 ul 

34 


In the tables for peak correlation parameters, entries ■ 
above the diagonal are for M(u,v), M(u, v) is given belov/ the 
diagonal and M(u) is along it. 

From the data available, it is conjectured that minimum 

odd cross- correlation which can be achieved with M-length 

c 

p-NRQFs is ^ i tends to remain constant as N 

increases unlike the case of m-sequences, where this ratio 
asyirptotically goes to zero amongst certain sets of sequences. 
Hence, periodic and aperiodic correlation performance is much 
worse vis-a-vis m-sequences. More data needs to be analysed 
to draw any further conclusions. . 

It is to be noted that NRQFs are not the only binary 
Boolean functions with uniform-magnitude sequency- spectrum. 
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Example 4 , 22 ; The repetitive quadratic form 

f ^ (x) = ® x^x^ ® ^2^3 another degree 

three function, 

f 2 (x) = x^X2X3 ® x^x^^ © X2X5 ® X3Xg 
are two non-NRQFs with uniform - magnitude sequency spectrum. 

The general class of these functions, called bent 
functions, is discussed in the next chapter. 



CHAPTER 5 


BENT FUNCTIONS 


This chapter/ the results of which were already avail vable, 
is included to complete the notion of NRQFs introduced in the 
previous chapter in which the work done had been completely 
independent of this. 

The notion of NRQFs 'is generalised to that of bent functions 
in this chapter. Sections 5.1 and 5.2 deal with their defini- 
tion and some of the properties. Section 5.3 presents the 
problems encountered in the enumeration of bent functions of 
a given number of variables. Section 5.4 goes on to the sub- : 
class of quadratic bent functions/ their graphical representa- 
tion and their exhaustive enumera'tion. Auto- correlation of bent 
functions is derived in Section 5.5. 


5.1 , DEFINITION [MacWilliams and Sloane/ 1977; Rothaus/ 1976^ 

Let f (x) be a binary Boolean function of n =2s variables. 


Let ^f (x) = l-2f(x) i.e. / f (x) in 
and F = H ^f be the u-WHT of '^f. 



Then f (x) is called a bent function iff each WHT coefficient, 
F^/ is +2^^^ (+.2®). 



112 


The already encountered NRQFs satisfy this condition. 
An example of a noh-NRQF bent function is: 

Example 5. 1 : Let f(x) = x^x 
function of 4 variables. f, and F are as in Table 5,1. 

Table 5.1 


WHT of f(x) - x.x.®x„x^@x^x^ 

1 1 2 3 3 hr 


X 

( x ^ 


» ^ 

2 ^ 1 ^ 

f 

Cf 

F 

0 

0 

0 

0 

0 

0 

1 

4 

1 

0 

0 

0 

1 

0 

1 

4 

2 

0 

0 

1 

0 

0 

1 

4 

3 

0 

0 

1 

1 

1 

-1 

,/! 

— 

4 

0 

1 

0 

0 

0 

1 

4 

5 

0 

1 

0 

1 

0 

1 

4 

6 

0 

1 

1 

0 

1 

-1 

4 

7 

0 

1 

1 

1 

0 

1 

-4 

8 

1 

0 

0 

0 

0 

1 

4 

9 

1 

0 

0 

1 

0 

1 

-4 

10 

1 

0 

1 

0 

0 

1 

4 

11 

1 

0 

1 

1 

1 

-1 

4 

12 

1 

1 

0 

0 

1 

-1 

-4 

13 

1 

1 

0 

1 

1 

-1 


14 

1 

1 

1 

0 

0 

1 

-4 

15 

1 

1 

1 

1 

1 

-1 

_4 


^ 2^3 


x^x^^ be 


{o,.} 


valued 


It is easy to see that this bent function, called 'a 
quadratic form, is an Affine Transformation of the p-NRQF, 

g(x) = 3C^X2 ® ^3^4 

(Two n variable binary functions f(x) and g(x) are said to be 
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related through an affine transformation iff 

f (x) = g (x A ® B) , where A is an n x n invertible 
binary matrix and B is a 1 x n binary vector) . 


In the above case. 



B (0 0 0 0) 


g(x A ® b) = (x^ ® x^) X 2 ® ^ 3^4 

= X-X„ © x^x_ © X-,X, 

1 2 2 3 34 

= f (x) . 


5.2 PROPERTIES OF BENT FUNCTIONS 

Let f (x) be a bent function of n variables, 

( 1 ) n is even' i.e., bent functions exist only for an even 

number of variables, n = 2 s and if n^ 2 , then deg H n. 

( 2 ) 2~^ times the u-WHT of a bent function of n =2s variables 

is also bent. 

F = H '^f is u-'WHT of f. 

F. = + 2® , i = 0, 1,..../ 2 ’^-!. 

X 

So, = 4; 1 = (- 1 ) , where f^(i) is another Boolean 

function. 

So, 2 “® F defines another function f^. 





Let = 1 - 2 f ^ 

NOW/ if = H / 

then/ F^(i) = +2®. 

Hence/ is also a bent function. 

f (x) has a natural partner in through its WHT. It is 

possible for certain bent functions that f and f^ are same. 

i.e./ H ^ f = 2® ^f (5. l) 

It is obvious that such functions are the eigenvectors of the 
Hadamard matrix H. 

Conjecture; All the p-NRQFs of ru;, variables are eigenvectors 
of the N X N (N = 2^) Hadamard matrix H. 

This has been verified for 2/4/6 and 8 variable p-NRQFs 

Exciirtple 5.2 : f, P and f^ for the three p-NRQFs of 4 variables: 

I 

f(x) = ® ^3^4 

^f(x) = x^x^ ® ^2^4 

3 

f(x) = ® ^2^3 shown in Table 5.2. 

. (3) F = H 

Fo = ± 2®. 

Let Z = Number of zeros of f 
= Number of ones of f. 
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Tgble 5.2 

f , F and for (x) , (x) and (x) : 

(x) = X ® x,x^ 

^f(x) = X^X 3 ® X 2 X^ 

3 

f (x) = X. X, ® X„X- 

14 Z O 


X 

(X4X, 

r^: 

X ) * 

21^ j 

- 1 


1 

F 

-d ; 

. - 1. . . 

f 


2 ' 

f « 
-d 1 

I. 



3 . 

^d 

0 

0 

0 

0 

0 

0 

4 

0 

0 

/I 

0 

0 

4 

0 

1 

0 

0 

0 

1 

0 

4 

0 

0 

4 

0 

0 

4 

0 

2 

0 

0 

1 

0 

0 

4 

0 

0 

4 

0 

0 

4 

0 

3 

0 

0 

1 

1 

1 

-4 

1 

0 

4 

0 

0, 

4 

0 

4 

0 

1 

0 

0 

0 

/I 

“x 

0 

0 

4 

0 

0 

4 

0 

5 

0 

1 

0 

1 

0 

4 

0 

1 

-4 

1 

0 

4 

0 

6 

0 

1 

1 

0 

0 

4 

0 

0 

4 

0 

1 

-.4 

1 

7 

0 

1 

1 

1 

1 

-4 

1 

1 

_4 

1 

1 

-4 

1 

8 

1 

0 

0 

0 

0 

4 

0 

0 

4 

0 

0 

4 

0 

9 

1 

0 

0 

1 

0 

4 

0 

0 

4 

0 

1 

-4 

1 

10 

1 

0 

1 

0 

0 

4 

0 

1 


1 

0 

4 

0 

11 

1 

0 

1 

1 

1 

_4 

1 

1 

“• ‘"X 

1 

1 

_4 

1 

12 

1 

1 

0 

0 

1 

_4 

1 

0 

4 

0 

0 

4 

0 

13 

1 

1 

0 

1 

1 

-4 

1 

1 

-4 

1 

1 

-4 

1 

14 

1 

1 

1 

0 

1 

-4 

1 

1 

—4 

1 

1 

_4 

1 

15 

1 

1 

1 

1 

0 

4 

0 

0 

4 

0 

0 

4 

0 


So, Zg - = + 2 and 


Zo + Z^ = 2 


2s 


r 22 s~1 ^ 2®'^ 


Thus/Z^ 


Zi H 


o2s-l nS-1 

2 f - 2 


2^®“^ - 2®“^ 


22s-1 ^ 2®”^ 
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Hence, a necessary condition .for f(x) to be a bent function 
is that its number of ones and zeros should be one of the above 
sets of values. 


(4) Any affine transformation of f (x) 
function. 

2-1 <m, x> V 

P ROOF t F(m) = S (-l) (- 1 )^'^’^'' 

1=0 

= +2^ 


results in another bent 


Let g(x) = f(x A ® b), be an affine transformation 


of f (x) . 


So, 


So, 


(m) 


N-1 

s 

i=0 




Let k = i A e B 
or i = (k ® B) A~^ 


G(m) = 2 


N-1 

r 

k=0 


(k®B) N"> (-1) f (k) 


<m',k> <m, BA“^> , ^^fCk) 

S (-1) •(-!) — 

k _ ^ T 

where m' =.-m(A ) 


(± 1). s (-!)<'"' , 


as 


(+1). 


(+ l).F(m'), as m varies from .0 to N-1, m' also 
takes all these values. 
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Thus, G(m) = + 2 ® . 

Hence, g(x) is also a bent function. ^ 

(5) Complement of a bent function is also bent. 

Properties 4 and 5 can be used to find bent functions of 
a given number of variables. 

5.3 NOTE ON THE ENUMERATION OF BENT FUNCTIONS 

Unfortunately, there is no kno•^^7n elegant method of enumera- 
ting all bent functions of a given number of variables. Property 
4 holds that affine transformations divide the set of bent 
functions into various equivalence classes. But, it does not 
give the representative function of each class whose affine 
transformations are the rest* 

In general, to find this representative function, a brute 
force approach using the definition of bent functions can be 
adopted. 2'^-point binary functions vjith weights (2^'® ^-2“^ 
are selected and their u-WHTs found one-by-one. If a bent 
function is encountered, all its affine transformations are 
grouped together and deleted from the set. Complements of these 
functions are also bent. 

Number of affine transformations for n = 2s variables 

= 2"^. TT^ ( 2 ""- 2^), and 
i=0 

2 jL s 1 2 ^ ' 

number of binary functions with weight (2 -2 ) = ( 
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Hence/ this method becomes impractical as n increases 
because number of possible functions to be checked increases 
e 3 <ponenti ally. 

If generalised bent functions are defined as real valued 
2^-point(N =2^/ n = 2s) functions with uniform-magnitude 
sequency spectrum/ then these can be exhaustively enurrierated. 
If f is such a function. 


r . . ,1 


H f=(N)"^ £ / v/here c[ +l,-lj - valued function. 


■F = 


(N) 


'2 = 


H g 


(5.2) 


2 ^' f 1 _ n 

As f is real-valued/ all possible 2 ' "5 4-1/ -Ir valued 2 -point 
functions £ will result in distinct generalised bent functions 
f of n variables. This set will be complete. 


Example 5.3 ; Consider the generalised bent functions of 

2 ^ 

2 variables. Number of functions a is 2 i.e. 16. The 
16 bent functions computed using (5.2) are: 


So 

1 

1 

1 

1 


2 

0 

0 

0 


1 

1 

1 

-1 


1 

1 

1 

-1 

22 

1 

1 

-1 

1 

^2 

1 

-1 

1 

1 

23 

1 

1 

-1 

-1 

23 

0 

0 

2 

0 

24 

1 

-1 

1 

1 

£4 

1 

1 

-1 

1 

25 

1 

-1 

1 

-1 

-5 

0 

2 

0 

0 

26 

1 

-1 

-1 

1 


0 

0 

0 

2 

27 

1 

-1 

-1 

-1 

-1 

-1 

1 

1 

1 



119 


3 -1 

1 

1 

1 

4 

1 

-1 

-1 

-1 

Sg -1 

1 

1 

-1 

fg 

0 

0 

0 

-2 

2l0 

1 

-1 

1 

”10 

0 

-2 

0 

0 

2U 

1 

-1 

-1 

~11 

-1 

-1 

1 

-1 

2i2 -1 

~1 

1 

1 . , 

h" -12 

0 

0 

-2 

0 

2i3 

-1 

1 

-1 

- ^13 

-i: 

1' 

-1 

-1 

3i4 

-1 

-1 

1 


-1 

-1 

~1 

1 

%5 

-1 

-1 

-1 

45 

-2 

0 

0 

0 


Again, for these computations, oniy . representative functions 
of each equivalence class need be computed and the rest are 
their affine transformations or complements* --For instance, 
computation of f^, f^ and fg is ..sufficient ..for finding 

out the rest. . 

Hence, in general, the enumeration of bent functions is 
still an open problem* However, a class of bent fxanctions, 
called quadratic bent functions, is anumerable. 

5,4 QUADRATIC BENT FUNCTIONS ’ ' ' 

Bent functions consisting of ;quadratic (degree' 2) terms 
only are called quadratic bent functions (QBFs). It is possible 
to find the total number of QBFs of a given number of variables 
in general. 

Any quadratic form (not necessarily bent), f(v), of n 
variables can be ejcpressed in terms of an n x n upper triangular 
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binary matrix Q. q. . is one iff v.v. is a term of f(v). 

^ J J 

Example 5.4; . Consider the QBE, f (v) , of 4 variables, 
f(v) = 


Alternatively, it is e 3 <pressible as; 



0 10 0 


^1 

r 


0 0 0 0 



f (v)=: 1 Vt V„V-,V. 





1 12 3 4 


0 0 0 1 


^3 


; 0 0 0 0 


J^4, 


T 

= y Q y , 


where Q (v) = 

Similarly, 

S 


0 

0 

0 

0 

0 

0 

0 

0 


10 0 
0 0 0 
0 0 1 
0 0 0 

10 1 
0 10 
0 0 1 
0 0 0 


represents 


g (v) = V Q y"^ 


= ViVj » ® V3V^ 


g(v) is not a QBE. 

Further, a quadratic form, f(v), can be characterised by 
a symmetric matrix, 

S = Q ® Q^' 

S generates another form which is symmetric bilinear, called a 



symplectic form/ the matrix S being termed as symplectic 
matrix. This form, in general, for a given S is : 


S (U/ v) 

= U S 

y 

T 



r- -1 

i u . . . u 


— 





11 nj, 





V 


["l--- "n] 





S 


(s. .) 
ij n xn 

,e 5.5: 

Matrix 

S 

for the 

QBE of previous example is: 



r 



-] 




0 

1 

0 

0 


S 


1 

0 

0 

0 




0 

0 

0 

1 




0 

0 

1 

0 

1- 

! 



r 



n 

i r 1 


S(u,v) =j 


0 10 0 

10 0 0 

0 0 0 1 
0 0 10 - 


V, 


V. 


V- 


V, 


U V 
2 1 


^l''2 






Evidently, each quadratic form has a unique Q and S 
representation and vice-versa. Whether a form is bent depends 
on the rank of S. 


Lemma 5.1; 

Let N(m, r) be the number of m xm synplectic matrices of 
rank r over GF(2). Then, 

N (m, r = 2h+l) = 0 and 

(2^’’-l) .. . (2^-1) 


N (m, r = 2h) 
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Lemma. 5.2 ; 

A cjuaciratic form, f(v), of m variables is a bent function 
iff its symplectic matrix S has rank r«m. i.e. 2h =m. 


Example 5.6 ; 

(i) Let f (v) = v^v^ ® "^2^3 ® ® ^ cjuadratic form 

of 4 variables: 

0 110 
10 10 
i ^ 110 1 

0 0 10 

Rank ( S ) =4. So, f(v) is a bent function. 


(ii) Let g(v) 

~ '^ l'^2 


0 

1 

0 

1 


V2V3 © v^v^ 


10 1 
0 10 
10 1 
0 10 


^3^4 


g(v) is not a bent function because rank ( § ) =2. 

For small values of m (say m ^ 6) , whether a quadratic 
form is a bent function can be found out by inspection from 
its S matrix. 

Lemmas 5.1 and 5.2 jointly give the total number of 
possible quadratic bent functions of a given number of 
variables. Number of quadratic bent functions of m = 2s 
variables is the niimber of binary m x m symplectic matrices 
with rank r = m which is: 
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For m = 4 , N = 28 

For m = 6 , N = 1388S 

For m = 8 , N = 112881664 

These do not include the functions with complementation of 
variables. 


Lemma 5.3 ; 

Any quadratic bent function of m variables is an affine 
transformation of the principal form 

fp(v) ® '^3'^4 ® '''’’rn-i ''^m ® ' ee^|0,l| 

e.g. bent functions of m=4 variables are all quadratic forms. 
They can be divided into four different classes having two, 
three, four and six quadratic terms respectively, each being 
expressed as an affine transformation of the principal form. 

<l) f(v) =s V V ® v„v,, (3 typesi 

J.. <iL w 

o 

( 2 ) f (v) = y^V 2 ® V 2 V 3 ® '^ 3'^4 (12 types) 



^f (v) = (A V ® B ) . 
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3 

(3) f (v) = ® ® '^3’^4 ® types) 








1 

0 

0 

0 


0 

0 

1 

1 

0 

B = 

0 

0 

0 

1 

0 


0 

0 

1 

0 

1 


0 


(v) = (A V ® B ) 

4 

(4) f(v) = v^V 2 ®v^V 2 ®v^v^^®V 2 V 2 ©V 2 V^®v^v^ (l type) 


“ 



10 0 0 


1 

110 0 


0 

0 110 

B = 

1 

0 1 0 1 


1 


^■f (v) = (A y © B) with e = 1 

The graphical representations of these 28 bent functions are 
given in Table 5.3. 


5.5 AUTO- CORRELATION PROPERTY OF BENT FUNCTIONS 
Theorem ; A Boolean function f(v) of m=2s variables is a bent 
function iff its dyadic auto- correlation is of the form: 

r 2^® , t =0 

b.Ct) =< 2s 

0 , t = 1, 2, . . . ,2 -1 . 

'-4. 


PROOF: 


(i) Let f(v) be a bent function of m = 2s variables 
F = H ^ f 

* O ' ^ 

So# F(x) 2 / 1 = 0 / l/*** 2 ' —I* 


NOW bj = 2^-\h ( ^ f . ^ 


E ) 



Ta-i>*e 


(X)( 3 



(E)(/atyfe^P 




(12} (1 type} 





12^ 


So, = ~~ 

' -f ^22 


= H. [: 


Thus, b.(t) =>1 2 , t-0 

lo , t=l, 2, . 


.2s-l 


(ii) Let f(v) be a Boolean function of m = 2s variables v;ith 
dyadic auto- correlation of the type given above. 


so, [_2^® 0 q] 


1 2 


Or, = 2^® b. 


M iif 


[ 22 ".... 22 § 


SO/ 08.clri Fj^ “ i 2 / 


Hence, f(v) is a bent function. 


Q • F • D • 



CHAPTER 6 


CONCLUSION 

We began with an integrated approach to correlation, 
linear systems and the system-specfic transforms. Design of 
signals with good correlation properties is posed as an 
important problem within this overall perspective. Sequences 
with ideal and asymptotic ally-ideal cyclic correlation proper- 
ties have been reviewed from the linear cyclic- shift invariant 
(LCSIV) systems' and discrete Fourier transform (DFT) points 
of view, thus, establishing necessary and sufficient conditions 
for any sequences to possess these properties. In fact, these 
results esqjound criteria which aid the synthesis of generalised 
sequences with desired correlation properties. However, binary 
sequences were dwelled upon more because of their obvious 
advantages in terms of generation and other properties. 

We passed on, further, to an exposition of the inter- 
relations between linear dyadic-shift invariant (LDSIV) systems, 
dyadic correlation and the Walsh-Hadamard transform, thus, 
preparing the groundwork for a study of sequences with good 
dyadic correlation properties. Non- repetitive quadratic forms 
(NRQFs), a class of binary Boolean functions, was dealt with 
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extensively as a set of binary sequences with ideal and asyitpto- 
tically ideal auto- and cross-correlation properties# respec- 
tively. The uniformr-magnitude Walsh-Hadamard transform of NRQFs 
was arrived at in a general way, and applied to their generalise 
dyadic correlation analysis. Ihis approach lead to some inte- 
resting results. In order to assess the periodic and aperiodic 
correlation performance of NRQFs# periodic and aperiodic corre- 
lation parameters were cortputed for some specific cases and 
compared with those of m- sequences. 

We also presented an integrated brief overview of bent 
functions# a larger class (than NRQFs) of binary functions 
with uniform-magnitude sequency- spectrum and hence# the ideal 
dyadic auto- correlation property. The vini form-magnitude spec- 
tral property suggests a method to obtain generalised bent 
functions. These were also presented, # The generalised 

bent functions constitute a complete set with ideal dyadic auto- 
correlation, This important observation marked the end of 
dyadic correlation analysis attenpted in this work. 


6,1 SCOPE FOR FURTHER WORK J 

Two problems related to the present work which have been I 

f 

investigated summarily are dealt with presently. First is the i 

design of binary sequences with specified dyadic auto-correla- > 

tions using the method of dyadic difference sets. This is analo-j 
gous to the synthesis of sequences with specified cyclic corx“e- 
lations using ..integer difference sets. The problem is i 
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introduced here, a detailed analysis being beyond the scope 
of this work. Section 6,1.1 introduces the method for cyclic 
sequences which is extended to that for the dyadic case in 
Section 6.1.2. 

Another problem of interest is the study of dyadic corre^ 
lation properties of 2-D arrays. It is introduced in Sec. 6. 1,3. 
Sections 6, 1,4 and 6.1.5 present two methods of synthesis of 
2-D arrays from 1-D binary sequences. Some interesting results- 
are derived with NRQFs and Walsh fxonctions as the 1-D coirponent 
arrays.' 


6. 1. 1 Integer Differences Sets and Periodic Correla-tion 

Let X = (x^x^. . be N-length sequence with 

k ones. The oub-of-phase periodic auto-ccrrela-tion is: 

N-1 X X. 

A(t) = 53 (-1) ^(-1) / t=l, 2,...(N- 1), (i-t) mod N. 


Let A be a constant for all t as is the case with pseudo- 
random sequences. Define m to be the number of places# in 
which snd ^ ^i-.-t^ have a one in common i.e., 

^i~ ^i-t ~ ^ exactly m values of i, O'^i^ N-l 

and for all t, 1 

So, -the arrangement of ones and zeros of 


as follows: 
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m 


(k-m) (k-m) N-irv- 2 (k-m) 






1 1****1 0««*»0 
1«««*1 0««*#0 Oi 


,0 ■ 

,0 


Thus, A = m(-l) ^(-1) ^ H-aCk-m) (-1) + (N-2 k+ m) 

= m-2k+2m + N-2k + m 
= N-4 k+ 4in 

For a given N and k, A is constant if there are exactly m 
overlapping ones between and { ^i. for all t. 

This is possible if each non-zero cyclic shift, t, is 
expressible as a difference mod N between each of exactly m 
pairs of numbers mod N, each number being i such that Xj^ = 1. 

Example 6,1 ; Let x be a 7-length sequence: 


0 

12 3 4 

5 6 

1 

0 0 10 

1 1 

is to 

be noted 

that 

t 

Differences 

1 

(0-6), 

(6-5) 

2 

(0-5), 

(5-3) 

3 

(3-0), 

(6-3) 

4 

(3-6) , 

(0-3) 

5 

(5-0), 

(3-5) 

6 

(6-0), 

(5-6) 


Ntombers in each pair constitute the set‘^0, 3, 5, 6^, every 
element of which is an i such the x^ = 1 in x. 

The problem of synthesizing sequences with given' N,k and 
A is related to the existence of integer difference sets with 



parameters (N, k, ra) . 


Integer Difference Sets CMeetham, 1969J; 

A set of k distinct integers (d^ / d 
an integer N is called an integer difference set D with para- 
meters (N, K/ m), if every integer b ^0 (mod N) can be ej^sressec 
in exactly m ways in the form d^-d . 5b (mod N) , v/here d,>d. 
belong to D. 

Hence, if for desired A, k , and m, such a difference set 
is found/ then the seguence with out-of-phase auto— correlation 
A can be formed. 


2 / . . . / dj,) modulo 


Example 6. 2 i Consider a 7— length seguence, A is desired to 
be -1. 

SO/ -1 = N - 4k + 4m. 

Choosing k = 4 , m = 2 


There exists a difference set , D, mod 7 with parameters 
(7,4,2). If D = -{o, 3, 4, 5 } with mod 7, 


then l54-3e5 - 4 


2 =5 - 3 =0 - 5 

3 =3 -0=0-4 

4 =0 - 3 =4 -0 

5 =3 -5 =5 -0 

6 =3 -4 =4 -5 


mod 7 


Assigning a ‘1’ to position numbers in D and *0' 



i 


to the rest 


we get. 


0 1 2 3 4 5 6 

1001110 which has the desired correlation: 
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7 -1 -1 -1 -1 -1 -1 . 

Hence^ each such D with parameters (N, k, m) corresponds to a 
binary sequence with out-of -phase auto- correlation (N-4k-4m). 


1. 2 Dyadic Difference sets and Dyadic Correlation 

Let x = (Xq. . where N = 2^, be a binary -^0,1^ sequence, 

^ ' t =1^2, (N-1) 

i=0 


is the out-of-phase dyadic auto- correlation. Again, let k be 
the ones in x and m be the nvimber of places in which ^x^^and 
its dyadic shifted version have a one in common for 

all t, 1 ^ t;^ (N-1) . 

So, b(t) = N - 4k + 4m, 


For desired b, N and k, a difference set with mod 2 (®) as the 
operation is required such that each i^O (mod N) can be 
represented in exactly m ways as the ® sxim of pairs of its 
elements. 


Dyadic Difference Set ; 

A set of k distinct field elements {^1 \] over a 

finite field of N elements is called a dyadic difference set, 
D, with parameters (N, k, m) if every field element b^^O can 
be e: 5 ^ressed in exactly m ways in the form d^® fh, where 
dj^, dj belong to D. 


1 

!' 


i; 


!' 


[: 
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One possible way of choosing a finite field of N elements 
to represent position numbers from 0 to (N-l), with ® as the 
operation, is using primitive irreducible polynomial over GF(2) 

of degree n (N = 2^) . This generates an extension field GF(2^) 
with 9 ss th 0 fxsld opsxration* 


As an illustrative example, let N=2'^=16. 

Use the degree 4 primitive polynomial 
4 

2 - z ® 1 to generate the field. 

Let a be the primitives element. 

4 

SOjf oL S5 3. ® 2.* 

All the field elements, their chosen binary representation and 
the position numbers they represent are as follows: 


Field 

Binary 

Position 

Element 

Representation 

Niombers 

0 

0000 

0 

15 

X = a 

0001 

1 

a 

9 

0010 

2 

iCi 

a 

0100 

4 

o 

a 

1000 

8 

4 

a = a ® 1 

0011 

3 

5 2 _ 

a = a w a 

0110 

6 

6 3^2 

a = a a 

1100 

12 

7 3^ 

a = a ® a ® 1 

1011 

11 

8 2 - 

a - a ® 1 

0101 

5 

9 3 ^ 

a = a w a 

1010 

10 

XU A ^ 

a =0 w a © 1 

0111 

7 

11 3 ^ 2 _ 

a s= a w a ® a 

1110 

14 


, . contd 



15 

13 

9 


a 

a 

a 


12 

13 

14 


= a^ ® a^ ® a 0 1 
= a ® a ® 1 
= a^ ® 1 


1111 

1101 

1001 
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For ideal auto-correlation, b(t) = 0 and let k = 6 if 
N Si 16 (taking a cue from bent functions) . 

So, 0 = 16 - 4 (6) + 4 m 

or m s= 2, 

Thus, a dyadic difference set with six elements is to be chosen 
such that all the non-zero field elements are es^jressible as ® 
sum of pairs of its members in exactly two ways. 

One of the possible choices isi 

^ f 5 10 7 14 13 111 

D=<a,a ,a,a ,a ,a j 

t 

It can be verified that: 


1 

ss 


ss: 

5 a, 10 
a ffi a 

SS 

10 

a 

0 

5 

a 

2 

ss 

a 

5S 

al« * a’ 

=: 

7 

a ® c 

14 

i 

3 

sz 

4 

a 

ss; 


ss 

11 

a 

9 

a"-' 

4 

s 

2 

a 

ss 

a“ * a« 

ss 

13 

a 

9 

al4 

5 

ss 

a® 

ss 

7 11 

a ® a^ 

ss 

all 

9 

7 

a 

6 

ss 

5 

a 

ss 

a’ ® a« 

ss 

al3 

9 

7 

a 

7 

s 

alO 

ss 

a“ » 

ss 

14 

a 

9 

11 

a 

8 

r= 

3 

a 

ss 

a® ® a^^ 

ss 

all 

9 

5 

a 

9 

s: 

a^^ 

ss 

aW « a“ 

ss 

11 

a 

9 

alO 

10 

ss 

9 

a 

ss 

a“ « a« 

ss 

13 

a 

9 

10 

a 


. . . contd. 
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11 


7 

a 


= a^^ ® a^ 

12 

=■ 

6 

a 

= a’ 9 

= a^" ® a'7 

13 

ss 

13 

a 

= a^ ® a"^ 

= a"^ ® a^ 

14 


11 

a 

= a^° ® 

= a^^ ® al° 

15 


ai2 

= a^ ® a^^ 

= a^^ © a^ 


Assigning a ' 1 ' to the position numbers corresponding 
to the elements of D and a 'O' to the rest/ we get the segtience; 


n ^15 ^ ^4 2 8 5 10 ^3 ^14 ^9 1 J ^13 ^11 12 
ua aaaaaa aa aaaa a a 


0 0 


000011 01 010 1 


0 


This has 'O' out-of-phase dyadic auto- correlation as 
desired. Incidentally, this is also an NRQF : f(x) =^y:^yi^9yi^yiy 


This outline suggests that if dyadic difference sets 
with parameters (N, k, m) can be constructed, then sequences 
with desired dyadic auto-correlation can be designed. Using 
GF(2^) as the finite field, all the primitive degree n poly- 
nomials have to be tried. Another method outlined by 
(McFarland, 1973^ can also be tried. There can be many dyadic 
difference sets arising out of a given degree n primitive poly- 
nomial. This also suggests a possible method for enumeration 
of bent functions of n variables. Each dyadic difference set 
with CN, k, m) such that b(t) = 0, t=l, .,.(N-l) and 
k = 2^”^ - represents a binary sequence which is a bent 

function. In other words, there is a one-to-one correspondence 
between such difference sets and bent functions. 
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6. 1. 3 Dyadic Correla'blon Properties of two*Dimensional Arrays 
Let C(x, y) and dCx, y) be two N x M (N =2^, M = 2^) 
2-dimensional (2 -d) real-valued arrays. Dyadic correlation. 


(i,j), is e3?pressed as: 

N-l 

M-1 

Rcd (i,j) =. r 

c(x, y) D(x0i, y®J) 

x=0 

y=0 


i =0,1, (N-l), 
j =0,1,..., (M-1) 


If c(x, y) = dCx, y), 

N-l M-1 

R (i,j) = ^ ^ C(x, y) C(x ® i, y ^ j) 

x=0 y=b 


Nature of the correlation function, R(i,j), depends on 
the method used for construction of the 2-D arrays j^Calabro 
and wolf, 1967] , The obvious choice falls on 1-D arrys having 
good correlation properties as basis of extension to the 2-D 
case. Two methods of forming binary 2-D arrays are discussed 
with exartples using NRQFs and Walsh functions as building 
blocks. NRQFs are taken to be J +l,-l"}’ - valued. 


6,1,4 2-D Arrays from Dyadic Shifts of 1-D Arrays 

Let, f(u) and g(u) be two N-point (N =2^) 
valued functions. Consider two 2-D arrays, c(x*y) and D(x, y), 
constructed as follows: 

C (x, y) = f (x 0 y) , 

d(x, y) = g(x » y) , C and D are both N xN , 
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H.— , (i # J } 


N-1 N-1 

51 ir C(x,y) C(x ® i, y ® j) 
x=0 yi=0 


N~1 N-1 

= 51 ' 51 

x =0 Y =0 
N-1 N-1 

= 21 -X 

x =0 ^=0 


f(x®y)f(x@y®i •j) 
fCz) f(z ® i ® j) , z =x ® y 


= N b^(i ® j), where b^Ch) is the dyadic 

auto- correlation function of f(x) 
Rows of are dyadic shifts of b^(t), 

(i) If f (x) is chosen to be an NRQF, 
b£(t) = N , and 

Rp(i,j) is a diagonal matrix with this choice of f (x). 


(ii) If f (x) is chosen to be , the kth row of H (the N x N 
Hadamard matrix ) , then, 

= *> iik • 

So, R^(i,j) = 

.r2 

s= + N 
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N-1 N-1 

S'. S' c:{x,y) Dix 9 ±, y ® j) 

x=0 y=6 

N-1 N-1 

X' X f (x ® y) g(x ® y ® i ® j) 
x=0 y=0 

N-1 N-1 

S' S' f (z) g(z ® i ® j) , 2 =x © y 
x=0 2=0 

N (i © j) 


(i) Let f(x) and g(x) be NRQFs. Then, rows o£ R„- are dyadic 
shifts of Each row will have that many non-zero 

values as b^^ (t) . 


Example 6,4 : 



(t) 


Let f(x) = x^X 2 © and 

g (x) = x^x^ ® ^2^4 ' N = 2^ 

8 , t = 0, 6, 9 
J -8 , t = 15 

0 , otherwise 


16 


Dyadic shifts for which R^j^(i,j) is non-zero are 

given in Table 6.1, Figures in parentheses are the correlation 
values (normalised) . 

(il) Let f = and g ^ be the choice of f and 2* 

known that dyadic cross- correlation function between two 
Walsh functions is identically equal to zero. 


i.e. bgg(t) = 0, for all t. 

So, R^j^(i,j) = 0, for all i and j which is similar to 
the ideal cross-correlation for 1-D arrays. 
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Table 6«. l 
j* = i «t 


t 

i 

0(8) 


6(8) 


9(8) 

15 (-8) 

0 

0 


6 


9 

15 

1 

1 


7 


8 

14 

2 

2 


4 


11 

13 

3 

3 


5 


10 

12 

4 

4 


2 


13 

11 

5 

5 


3 


12 

10 

6 

6 


0 


15 

9 

7 

7 


1 


14 

8 

8 

8 


14 


1 

7 

9 

9 


15 


0 

6 

10 

10 


12 


3 

5 

11 

11 


13 


2 

4 

12 

12 


10 


5 

3 

13 

13 


11 


4 

2 

14 

14 


8 


7 

1 

15 

15 


9 


6 

0 

bnsider f ■ 

= and g (x) , an 

NRQF, as the next > 

bjg(t) = 

+ (N)^ 






fence, 

(i/j) * 

+ 

(N) 3/2 





s: 

i 

(23=) , 

N = 

= 2^® 


2-D Arrays 

as Termwise Product 

of 1-D Arrays 


2-D arrays can be constructed frcsn term-by-tem product 
of 1-D arrays. Consider, A(x,y), an N x M <N = 2"^, M = 2'”) 
array formed using an N-length array C(x) and an M-length 
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one, B (y) such that: 

A(x,y) = C(x) D(y) 


N- 1 M- 1 

Then, R (i,j) = ^ ^ A(x,y) A(x ® i, y ® j) 

x=0 y=0, 

N-l N-1 

= 51 21 C(x) D(y) C(x ® i) D(y ® j) 

x=0 y=6 

= (i ) b^ ( j ) , 


is the termwise product of auto— correlations of C(x) 
and D (y) . 


(i) If C(x) and D(u) are N- and M-point NRQPs, 

b^(i)=NS^0' 

~ ^ ^ jO 


Hence, R. (i, j ) = NM S/ . . ... 

A (x,j)/(0,0) i.e. It is 


zero for all 


non-zero shifts. This is similar to the ideal auto- correlation 


property of 1-D arrays. 


(ii) If C(x) and D(y) are N- and H-point Walsh functions, 

R^(i,j) = (+N) . (+M) for all i,j- 

= + NM. 

(iii) C(x), an N-point NRQF, and D(x), an M-point Walsh function, 

is another choice. . 

Then, R^(i,j) = N Sj_q. (± 

= ±™6i0 
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